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ABSTRACT 


The controls analysis based on a truncated finite element model 
of the 122m. Hoop/Column Antenna System focuses on an analysis of the 
controllability as well as the synthesis of control laws. Graph theore- 
tic techniques are employed to consider controllability for different 
combinations of number and locations of actuators. Control law synthesis 
is based on an application of the linear regulator theory as well as 
pole placement techniques. Placement of an actuator on the hoop can 
result in a noticeable improvement in the transient characteristics. 

The problem of orientation and shape control of an orbiting flexible 
beam, previously examined, is now extended to include the influence of 
solar radiation evironmental forces. For extremely flexible thin 
structures modification of control laws may be required and techniques 
for accomplishing this are explained. Effects of environmental torques 
are also included in previously developed models of orbiting flexible 
thin platforms, and a preliminary analysis of related thermally induced 
deflections and torques on thin structures is provided. Finally, it 
is concluded, based on an evaluation of the coupling coefficients in 
the rigid modal equations and the generic modal equations for the first 
seven modes of the Hoop/Column System, that when the system is performing 
within nominal specifications (surface deflections of the order of mm) , 
that the relative magnitude of the largest of these coefficients is at 
least one order of magnitude smaller than the principal coefficients. 
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I. INTRODUCTION 


The present grant extends the research effort Initiated in previous 
grant years (May 1977 May 1982) and reported in Refs. 1-8*, Techni- 
ques for controlling both the shape and orientation of very large in- 
herently flexible proposed future spacecraft systems are being studied. 
Possible applications of such large structures in orbit Include; large 
scale multi-beam communications systems; earth observation and resource 
sensing systems ;orbitally based electronic mall transmission; and as 
orbital platforms for the collection of solai energy and transmission 
(via microwave) to earth based receivers. 

This report is subdivided into seven chapters. Chapter II is 
based primarily on a paper to be presented at the 34th International 
Astronautical Congress and presents preliminary results on the controlla- 
bility and control law synthesis for the Hoop/Column orbiting large 
flexible antenna system. Graph theoretic techniques previously intro- 
duced in Ref. 7 and further described in Ref. 8 are employed to consider 
system concrollability for different proposed actuator arrangements 
(number and location). Control laws are then designed based on the 
ORACLS computer algorithm , primarily using linear quadratic Gaussian 
techniques. The mode shapes (eigenvectors) and frequencies of the system 
were obtained previously by finite element techniques and supplied by 
NASA-LRC,^° 

At the operational altitudes of the future missions involving 
large space structures, the principal environmental disturbance is 

*For references cited in. this report, please see list of references 
at the end of each chapter. 
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that due to solar radiation pressure. The effect of solar radiation 
(pressure) disturbance on a flexible uncontrolled orbiting free-free 
beam was examined in Ref. 8, and to the authors* knowledge represented 
a first attempt to Include such dlsf ,::bances in the system dynamics of 
a flexible structure in orbit. Chapter III of this report, based on 
a paper recently presented at the 1983 AAS/AIAA Astrodynamlcs Conference, 
extends this work to consider the closed loop dynamics of an orbiting 
flexible beam under the influence of solar pressure induced moments. 
Actuator control laws previously designed for the case where environmental 
effects were neglected are now re-evaluated especially for extremely flexi- 
ble beams. Shape and orientation control (about two different nominal 
orientations - local horizontal and local vertical) are considered. 

In Chapter IV, the effect of solar radiation pressure on the open- 
loop dynamics of a more complex structure, a thin homogeneous flexible 
square plate is analyzed. The mode shapes and frequencies of the plate 
are obtained numerically using a finite element computer algorithm. As 
in the case of the beam only small transverse deformations ate assumed, 

Solar pressure torques can also result on large flexible space 
structures due to the interaction of the Incident solar radiation on 
the thermally deflected bodies. Chapter V presents a preliminary formu- 
lation of the solar radiation heating effects on thin rectangular plates 
and beams. The amount of thermal deflection depends on the thermal 
gradient across the structure, the thermal conductivity of the material, 
the thermal expansion coefficient, and the thickness. Expressions for 
the thermal gradients can be obtained by simultaneously considering 
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the thermal radiation and heat conduction properties. These expressions 
can be used to develop equations for the thermal deflections and sub- 
sequently a model for the Induced moments on the (thermally) deflected 
structure. These moments will have to be considered In addition to 
those generated by the interaction of solar pressure on the freely 
vibrating system. 

Chapter VI describes the development of a computer algorithm to 
evaluate the relative magnitude of coupling coefficients In both the 
rigid and elastic modal equations of motion. These coefficients describe 
coupling between the rigid and flexible modes and also Intramodal 
coupling and are usually neglected when a finite element analysis is 
employed. The evaluation of the relative magnitude of these coefficients 
Is based on the numerical parameters for the Hoop/Column system. 

Chapter VII describes the main general conclusions together with 
future recommendations. The effort described in Chapters II, IV, and 
V is being continued during the 1983-84 grant period in accordance with 
our proposal submitted Jan. 1983.^^ 
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II. ON THE CONTROLLALILITY AND CONTROL LAW DESIGN FOR 
AN ORBITING LARGE FLEXIBLE AITTENNA SYSTEM 


Abicrteg 

Tb« cooerollAbillg^ ud gontrol l«v aTtjchcsxs 
b**«d on a Half aliusciic aod«l of eh« Boop/Coluaa 
orbiclng aacaua ayscaa is >*.3nsid«:si! . Cgaph chso- 
rsele esehniquss arc aaploysd eo analyza ehs coo- 
erollabilicy *>f cha systaa for posslbla pcoposad 
accuator arranfasnanu wbiob ineluda eorquara and 
pr^e aeruaeors aloog eba ease and a aiagla aerua- 
Co plaead on tba hoop aaaaebly. Oaea eoacroU- 
abliiey la astabliahad for a glvan eoebinacloa of 
oumbar aod locaelen of aeeuacora, sha synsbasla of 
coacrol laws la baaad on aa appLieaclon of tba 
linear reguliicor eheory aad also pola placeoase 
Cacbnlquaa. Za ganaraX. avrfaca torsloa and faad 
mast torsion are among sodas having cha longase 
cisa conataacs. Syscan traaalaoe parformacca is 
BoCleeably degraded whan the hoop sctuacor is noc 
laeludad. 

laeroduccioo 


Large flaxlbla orbiting syscaes have been pco> 
posed for posslbla use in eormmlcatiogs, elec- 
tronic orbital based oail ay stems, da cac cion of 
earth resources, and in solar energy coUaetion. 

The site and low walghe to area ratio of such sys- 
Ceaa dictate that sysrca flaxibilicy is now cha 
amio eonaideration in cha dynamics aad control prob- 
lem as eontr-i •. with She inherently rigid nacura 
of earlier . ■ -.aft lystams. For such large 

flcr' >1^ .CCS both orientation and shape con- 
trol vUJ. ;*ten be required. It is the purpose of 
Che presenc paper to consider the concroHabilicy 
and the eubsequent design of control laws based on 
e model of a proposed orbiting ancitsna system - 
Che Boop/Coluan aystaa. 

The Hoop/Coluan antasne system is one of the 
configurations under consideration for use in the 
foCure Bulcl-beaa Land hobila Satellite System^, 
designed to provide point eo point eocBunicatioas 
for Z20,000 subscribers across the V.S. in the 
aid 1990 's. The system is based on a large gso- 
synehreuous relay antasna and a n«abcr of eoblla, 
Ssrth-basad receivers. In ordar to achiava the 
required SF performance a pointing accuracy of 
(0.03 - 0.10} degree ?XS and a surface (antenna) 
accuracy of 12sa SMS will be reqxiired. The Boop/ 

* ColtE=a antenna ayscea^, depicted in rig. 1, in 
deployed configuration, ecntiins the deployeble 
(telescoping) nest system connected to the hoop by 
support cables under tcasinn. 




The hoop esntalns 43 rigid sections, to be deployed 
by nocor drive units. The daairad shape of the IT 
raflactiva naan is produced by « aeeecory d:r; ng 
surface using surface concrol cables, lha n .Isw- 
tiva aaah is eonnactad to tba boop by Juar: t 
graphite stringers. At one end of Cha nast c.ha 
alactronlc faad aisanbliss era positlsnad, wharaaa 
at the other end are tha principal solar arrays 
eonnactad to tha cain bus based control. 

Bacenely graph theoretic tachniquas^ were used 
Co study eontroUabilicy of linear syatans which 
could zaprasant large flexible orbiting systems 
with ifihartne damping. It was seen chat the tart, 
rente deficiency in the stiffness nacrix dictates 
Cha caq'sirsd r.fctbar as wall as tha location of cha 
aetuatora rsquirad for ccntroUsbility, wharess 
cha presence of the damping nacrix does not in- 
flusnea tha required sueber of actuators but pro- 
vides graatar flexibility to the availabla acruator 
locations for which the system is concrollsbla. 
Specific examples vers based on a previously de- 
veloped ncdel^ of an assumed homogeneous shallow 
Chin spherical shell In orbit, where both orlenca- 
clon and shape control were assumed to se provided 
by point actuators placed on the shell’s surface. 

It is the objective cf this paper to apply 
these techniques, first to analyze the controll- 
ability of Che Boop/Coluoa system based cn a possi- 
ble arrasgenant of actuators, aad, once controll- 
ability is assured, to consider various ways to 
develop control lavs. A cacenc related treatment^ 
of Che preliminary controls synthesis for this sys- 
tem was based cn a a/odel which contained three 
rigid-body rotational modes and a number of flaxu- 
I ral modes, but ignored the three rigid transla- 
clonal PLO^s. Actuators in the fora of torquers 
' only, were assened to be placed only alcsg t.4e nast 
assembly, aad ccntroUabili'ty was established in- 
directly through Duaerlcel neans. Control syn- 
thesis was acerrplisbed based on linear quadratic 
Gaussian techniques and also a related linear 
analysis assuming the collocation of actuators and 
sensors. 

II. ?!aiheeatical Model 

In the present paper, a finite elenant model of 
the Boop/Columo systen chat incltides all six rigid 
nodes and cha flaxiufal modes is used, “or the pur- 
poses of this study the systen is nodelled using 
112 node (grid points) with each node having a maxi* 
mun of six dagrees of freedom.® Thus, the Lin- 
earized dynamics of the Boop/Colustn antema systam 
can be described (in the absence of structural cr 
other damping) by; 

« -4 s - 7^ 
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fiS' lo 1‘be Boop /Colusa cnteima s^sceo. 
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iihar« 


m 



M Is ths 672x672 mass/ Inertia matrix 
K is ths 672x672 stiffness matrix 
X is ths 672x1 matrix consisting of the 
displscsotncs and rotations sc the nodal 
points 

F Is the control vector of dlasensiooslity 
672x1 


I After .<«lng the crsnsfacmacion between Che 
> nodal coordinates, q, and X given by 

X-<dq (2> 

i squacion (1) con be retirltCea as: 


q + 4>^EJq - 


(3> 


:'The left side of equation (3) con be expressed in 
Cerms of diagonal matrices, based on the pro- 
t parties of the eigenvalues and associated eigen>- 
vectors, according to: 


r “1,3 5 + c' Ki j ^ 


: where 




■ dlag. Ca^] - C'o^ ^ ] 

■ dlag. CX^] - C*Ki^ ] 

“ i^ genecollsed (modal) mess 
“ i^^ generalized stiffness. 


The first thirty^four (including the six rigid 
modes) eigenvectors end eigenvalues derived from 


the NASXBAX software for the system described by 
equations (1) and (2) were supplied;? however, 
only the first 26 modes had appreciable non-sero 
values 'or the corresponding modal aaases. 


The Influence of a finite number (m) of con- 
trol actuators or torquers is Incorporated into 
th* Bwdel by 

f' i ; 

: Fg - BO (j) ; 

i i 

jwhere 

' , B * 672»i control influence matrix 

0 * Bxl matrix associated with the control 
' vector 


jibe control influence matrix, B, (672xa) is formed 
las followe: 

' Zf there is on actuator chat influences the 
node (1< i < 112) in the Jth direction 
(1 < i < 6), then B(K,L)-1 where K - (l-l)x6+j, 
j and L " number of the actuator. Thus Che B matrix 
^ eooslscs of zeroe and ones, showing the influence 
i of ch« force octuecors on Che translational de- 
grees of freedom of the nodes, and also< that of 
ch« torque actuatore on the rotational degrees of 
freedom. 


The model considered here consists of six 
rigid body modes (3 translation -t- 3 rotation) and 
the first seven flexible modes. The values of the 
generalized masses and stiffnesses are listed in 
Table 1. 


Table 1 - Uoop/Coluen single layer surface 
nodal tlgenvalues 


Mode 

Fretuency 

Generalized 

Generalized 

No. 

B?. 

mass, ^ 

(Ib-see^/la) 

X, Stiffness 
(tb/ln> 

1 

0.0 

16.44388 

O.U 

2 

0.0 

8.923020 

0.0 

3 

0.0 

7.349353 

0.0 

4 

0.0 

9.704152 

0.0 

5 

0.0 

2.940652 

0.0 

6 

0.0 

8.418909 

0.0 

7 

.1188347 

153.1373 

85.38542 

8 

.2142455 

i. 232954 

9.482657 

9 

.2709358 

3.073094 

8.907021 

10 

.5063228 

.3046446 

3.083247 

11 

.7288725 

1.992988 

40.88663 

12 

.8897594 

723.5216 

22612.90 

13 

.9192313 

.6531203 

21.95405 


III. ControUabllitv Consideraticos 


Equation (4) con be rewrlctea with the aid of 
equaClon(S) as: 

q -<'a, ] C'K, ] q + C'a, ^^35 (6) 

Lh 


and then east into standard state space format with 
Che result: 



The number and location of the actuators ore 
determined after constderlng the concrollabllicy 
in addition to physical considerations. The dyna- 
mic system represented by equation (7) is con- 
trollable^ if and only if the pair 

-C'm, 3“^ C'K, ] and C'a, ]~^ <>‘^8 is controU- 

1 % Xh 

ble. 


The system matrix, S “ -C'a, 1 

for Che 13 mode model con be evaluated ax: 



( 8 ) 
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Th* digraph of Cho nacrix S Is shown In Tig. 2. 

For controllability^ cha actuators must be placed 
In such a way that they will lufluencn all the 
Biodes directly slnce> In this nodeli there is no 
coupling between the modes. The number of actua- 
tors will be dictated by the term rank (7) of the 
oatrl-x, 5. Because of Che six rigid modes (with 
zero eigenvalues)) the system matrix, S, has a 
rank deficiency of six; to augment Che term rank 
must have at least six linearly Independent 
non-zero columns, indicating a minimum of six pro- 
perly placed actuators la needed for controll- 
ability. 

A possible prsposed arrangement of a maximum 
of 13 actuators consisting of combinations of 
point actuators and a torquer, assumed to be 1^ 
cated along the mast and at selected points In the 
feed esaembly. Is deplctsd In Fig. 3. In addi- 
tion, e single point ectuaCor Is aesumed to be 
positioned on one of the rigid links of the hoop 
assembly. The actuators are assumed to be placed 
so that the vacloua modes Indicated In Table 2 
would be directly affected. 


Table 2 - Relationship between actuators 
and nodes directly Influenced 


Actuator no. Mode being affected 

(circled In Fig. 3) 


1 2 3 and 4 
3 

3 

7 

8 (Torquer) 

9 

10 

U 

12 

13 


Feed Mast Torsion (12) 
First Bending (about y 
axis) (B) 

First Bending (abduc x 
•xls) (9) 

Surface Torsion (10) 

Tsw (rotstlon about z axis) 
and First Torsion (7) 
Trsnslaclon along x and 
2nd mast bendieg (11) 
Translation alcng 7 and 
2nd mast bending (12) 
Translation along z 
Fitch (rotation about y 
axis) 

Roll (rotation about z 
axis) 


Possible sets of the minimum number of actua- 
tors Include (1,2,3,4.5,12) and (8,9,10.11.12,13). 
On Che contrary, the six actuators (1,2, 3, 4, 6. 7) 
are Insufficient to control the thirteen system 
nodes as seatea 14,15, and 26 in the digraph 
(Fig. 2) can not be reached from any of these six 
Inputs (under the assumption that any element In 
Che matrix which Is leas chan 10~^, and six 
or seven orders of magnitude less than the maximum 
entry, is created numerically equal so zero). 


tv. Central Law Svntheaia 

Onca controllablllcy has been established 
control law design Is based on: 1) application of 

the linear regulator problem from optimal control 
theory: and 2) pole allocation algorithms. The 
0RACL5° software package has been used to obtain 
Che control law gains and generate the required 
time histories of the actuator forces as well as 
the dynamic transient responses. 


X) Application of Che Linear Ragulacor Problem 
The control law of the fora: 
n - -FX (9) 


for the dynamical system 
X -AX+BO 


( 10 ) 


Is selected such that the .'oU.owlng performance 
index is minimized^ 


J - (x’^QX+U^RU) dt (11)* 

o 

where Q and R are positive seml-deflnlce and posi- 
tive definite weighting matrices, respectively. 


The gain matrix, F, Is given by: 

F “ R"Vp (12) 


where F corresponds to the positive definite solu- 
tion of the matrix Rlccatl differential equation: 

-PA-a’’p+PBR“ Vp-Q - dP/dt (13) 

The steady state solution (as c -* ■■) of the 
Rlccatl matrix equation can be obtained by solving 
equation (13) after setting dP/df-O. The steady 
state Rlccatl matrix solution results in a con- 
stant gain matrix which Is relatively easy to Im- 
plement as compared with Che time varying solution 
of the matrix Rlccatl differential equation. 

Tn order to guarantee conCrollabillty In Che 
event that one of the actuators might fall. It was 
decided to select t>ie minimum number of actuators 
la this study to be seven instead of the six pre- 
viously discussed. For seven actuators asstm»d 
(actuator numbers 7 to 13 la Table 2) a parametric 
study was performed shoving the effect of varying 
Q (1001 CO 10,0001) and R (I to 1001, where 1 Is 
the appropriately dimensioned Identity matrix) on 
the least dasq>ed mode of the closed-loop system 
(Fig. 4). It can be concluded that Q ■ 1000 I and 
R^l la a reasonable operating (starting) point in 
Che Q and R domain as the curve correspcr'''ing to 
R>I Is relatively less sensitive to changes In Q as 
compared with ocher values of R. 

As an example. Table 3 shows Che effect of 
varying the number o( actuators (from seven to 
thirteen) on the maxlmuo force amplitude required 
of any single aef’seor for the same set of Initial 
conditions In the modal coordinates and the s.'me 
state and control weighting matrices. As Indicated, 
the maximum peak force required of any actuator is 
reduced as the nui:d>er of actuators is Increased. 

The system transient performance Is Improved by in- 
creasing tho number of actuators (Table 4) . but a: 
Che expense of a alight increase In control effort. 
Ochar scuties showing the effects of: (1) removing 
three of the rigid modes from the model; or (2) 
ignoring Initial displacements In these nodes, 
could result in overly optimistic results in terns 
of control system requirements (T^^les 5 and 6) . 
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Fig. 2 Digraph of the aystem matrix. 



Fig. 3 Proposed atrangement of actoators Hoop/Coivra system. 

(uneireled numbers identify grid points from finite element analysis) 
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T«blB 3 - Sooaary of aiax'jinim Actuator forces(amplltudas) required 
Hoop/Coluan - varying cha no. of actuators (7-13) 

Initial Condition: q^(0)-0.01 1-1, 2.... 13, Q-IOOOI, R-I 


Max. Actuator 
Force Anplltudas 
Pounds 



0 

0 

fa 

Actuators 



7 

8 

9 

10 

11 

12 

13 

*1 







0.3500 

'2 




* 


0,1364 

0.0569 

'3 





0.3158 

0.3357 

0.3000 

'4 




0.3251 

0.2101 

0.1364 

0.0569 

'3 



1.303 

1.3029 

1.3027 

1.3027 

1.3027 

«6 


0.6632 

0.6631 

0.1831 

0.2898 

0.2847 

0.2833 

h 

1.694 

1.514 

1.315 

1.4734 

1.2798 

1.2781 

1.2761 

tg (in. -lb) 

0.109 

0.104 

0.104 

0.1004 

0.0124 

0.0U7 

0.0123 

*9 

1.1691 

1.1889 

0.287 

0.2276 

0.2369 

0.2369 

0.2369 

«10 

1.0602 

0.7380 

0.7380 

0.1673 

0.1601 

0.1573 

0.1573 


0.4209 

0.4089 

0.4326 

0.4199 

0.4055 

0.4061 

0.4061 

'12 

0.1208 

0.1206 

0.3521 

0.3521 

0. 3521 

0.3521 

0.3521 


0.2993 

0.1866 

0.1866 

0.0669 

0.0658 

0.0673 

0.0688 


ORIGINAL PAGE 19 
OF POOR QUALITY 


T«ble A - Comparison of closad loop poles - Hoop/Colusa with 7 and 13 actuators 


Q ■ 10001, R~Z 13 modes 


7 Actuators • 

Poles 

13 Actuators 

- Poles 

0-(Raal)-l/sec 

3b>- (loaginary) 

9- (Real) -1/ sec 

Ju- ( Imaginary) 

-0.0077 

5.5905 

-0.008 

5.5905 

-0.0077 

-5.5905 

-0.008 

-5.5905 

-0.0100 

0.7465 

-0.7045 

0.7788 

-0.0100 

-0.7465 

-0.7045 

-0.7788 

-0.4343 

0.4053 

-0.4179 

0.4544 

-0.4343 

-0.4053 

-0.4179 

-0.4544 

-0.8578 

1.7117 

-2.1916 

1.6704 

-0.8578 

-1.7117 

-2.1916 

-1.6704 

-0.9035 

2.7007 

-0.9806 

2.6637 

-0.9035 

-2.7007 

-0.9806 

-2.6637 

-3.0979 

1.666 

-3.6151 

1.7994 

-3.0979 

-1.666 

-3.6151 

-1.7994 

-0.4689 

1.3656 

-1.0058 

0 

-0.4689 

-1.3656 

-1.0117 

0 

-0.5583 

0.4892 

-1.0247 

0 

-0.5583 

-0.4892 

-1.0260 

0 

-1.0059 

0 

-1.0868 

0 

-1.0226 

0 

-1.0987 

0 

-1.0287 

0 

-1.2914 

0 

-1.2830 

0 

-1.5808 

0 

-1.5977 

0 

-3.9913 

0 

-2.4285 

0 

-5.4766 

0 

-3-42»jr 

0 

-6.2706 

0 

-6.9508 

0 

-10.019 

0 

-11.8231 

0 

-12.424 

0 

-18.1405 

0 

-30.186 

0 




V 


» • .s 


V. 


pr„. ., 


Table S - CoxparlaoQ of oaxlaum accuacor fcrce amplitudes 
Q • 10001, R-I, q^(0) - O.Ol, 


Max. Actuator 
Force Amplitudes 
Founds 

13 Actuators/ 13 Modes 

10 Actuators/: 

'l 

0.3300 


'2 

0.0S69 


'3 

0.3030 


'4 

O.OS69 ' 

0.2314 

*5 

1.3000 

1.3759 

'6 

0.2830 

0.2830 

*7 

1.2700 

0.3713 

fg (Inrlb) 

0.0124 

0.0171 

'9 

0.2360 

0.4031 

• '10 

0.1370 

0.2508 

^11 

0.4060 

0.1634 

'12 

0.3520 

0.2078 

'l3 

0.0688 

0.4115 


Table 6 ~ Comparison of maximum actuator force amplitudes 
Q - 10001, R-r 13 Actuators/13 Modes 


Max. Actuator q^(0)*O.01, i>l,2...I3 

Force Amplitudes 
' Founds 


qj^(0)-0, 1*1 CO 6 
qj(0)-0.01. 1-7 to 13 


fl 

0.3300 

0.3462 

S 

0.0569 

0.0324 


0.3030 

0.2923 

1 

0.0569 

0.0324 

* 

1.3000 

0.9846 


0.283 

0.3021 

^7 

1.2700 

0.8423 

fg (In.-lb) 

0.0124 

U.0124 

'9 

0.236 

0.1471 

* ^10 

0.157 

0.1566 

' hi 

0.406 

0.0327 

:j '12 

0.3320 

0.0976 

1 ^3 

0.0638 

0.C931 
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2) Application of the Pole Placement Algorithm 

Aa.an example of the application of the pole 
placement technlquot Fig. 5 llluatratea the tran- 
aient reaponae in key modal ampllcudea, as well 
aa required actuator force history for an Initial 
condition only In the first torsional mode. The 
pole placement algorithm^ requires that all of the 
poles be placed along a line parallel to the imag- 
inary axis. It is sesn that the requirement that 
each mode have a time constant of 100 sees, is a 
lest stringent requirement on Che control eyatem 
than for the cate depicted in Table 3 and 4 (13 
eceuatora) where the moat damped mode hes a time 
constant of approximately 1/30 sec. 


3) Effect of Removing Hoop Hounted Actuator 

If the seventh actuator Is removed from the 
system leaving a total of 12 working actuacora, 
thure la a deterioration in tha lowest damped 
modal Clma constants. The magnitude of the real 
part of the least damped mode Is reduced from 

O.OOB (Table 4) Co 0.0024 end there are now ten 
poles with amplitudes of real parts less thsn 1.0. 
The maxinuffl actuator forces required in response 
to Inlclel displacements of 0.01 in ell 13 modes 
ere not greatly different from those shown in 
Table 3 for the ease of 12 or 13 actuators. 

OD tha other band if a total of only seven 
actuators is now assumed (f^ 

there is considerable loss or perfonuace In Che 
least damped mods. The least damped modal time 
constant has now lncr.:ased to about 2000 sec as 
compared with 130 sec (as shown in Tsble 4} . It 
' should he noted chat the least damped modes 
I throughout this study correspond to both feed mast 
' toralon and surface torsion. Agein Che maximum 
eccuator forces requlrad in response to initial 
displacements of 0.01 in all 13 modes see essen" 
tially Che same order as shown In Table 3 for the 
i ease of seven actuators (including fy). 

i Additional results explained in Ref. 10, consi- 

i der the effect of varying the elements in the state 
penalty matrix so chat Che modes with poorest 
transient characteristics would be penalised more. 
There appears to be soma advantage in appropriate- 
ly aelecclng the "split” weighting elements so as 
to i apvove performance in the lowest damped modes 
without unduly increasing the actuator-force re- 
quirements. 


V. Concluding Comaents 

Generally, it can concluded that the use 
of the graph theory enables the designer to have 
an intuitive qualitative idea about where aecua- 
tora can be placed in order to achieve controll- 
, ability. In most cases simulated, surface torsion 
and feed mast torsion are among modes having the 
! longest time constants. Transient performance 
can be improved by appropriately changing the 
related elements In the state penalty matrices, 

, but usually at Che expense of control system ef- 
' ' I fort. The linear quadratic Gaussian technique 
offers more flexibility to the controls designer 
in attempting to meet the perfocmance requirements, 
^ ’ '^lle also maintaining propellant consumption at 
-'^...desired levels. System transient characteristics 
'* ^re degraded when the hoop-mounted actuator is 


not included. This degradation is more pronounced 
at Che number of (remaining) actuators is reduced. 
Tha control efforts depicted here for the transient 
times considered appear reasonable in terms of the 
over-all system mass and surface and pointing re- 
quirement. Before meaningful conclusions can be 
made regarding RMS pointing and surface accuracy 
requltements tha effects of both sensor and plant 
noise ei’id also environmental disturbances should be 
Incorpo' seed into the existing model and are sug- 
gested for further study. 
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POLE PLACEMEN! 

T - 100 sec. 

ntrriAL coNomoss; 



Time (secs) 

(All other nodes are excited very little.) 
TIME EESPOXSE - MODAL AMPLITUDES 


200 

TIME (secs) 


-H 

400 


600 


(Peak forces for ail ocher actuators are far less) 


'»odal amplitudes and 
required actuator force time history. 



II. B ADDITIONAL RESULTS ON CONTROL LAW DESIGN FOR HOOP-COLUMN USING 

LINEAR QUADRATIC CONTROL THEORY 

The concept of split weighting of different modes In the quadratic 
cost function Is used to further penalize modes 7 (combinations cf yaw 
and first torsion of the system), 10 (surface torsion), 12 (feed mast 
torsion) and their rates. These modes are penalized 10 times more 
than the rest of the modes and the resulting closed loop poles are 

tabulated in Table 7 together with the uniform weighting results. 

It Is observed that the time constant of the least damped mode (-0.0120) 
for the case of split weighting is improved without increasing the 
most damped mode (-30.1864) as compared to the uniform weighting case. 
The effect of split weighting Is also apparent in the maximum actuator 
force amplitude [Table 8], with the requirement of a slightly larger 
control effort. 

The actuator of the hoop (//7) is then removed; for the values of 

Q =» 1000 I and R = I the closed loop poles are compared with those for 

the case with the actuator on the hoop in Table 9. The corresponding 
required actuator forces are given in Table 10. The closed loop poles 
with the hoop actuators removed for a total of 7, 10, and 12 actuators, 
respectively, and with: Q ” 1000 I, R = I; Q = 10,000 I, R » I; and 
Q - 100, 000.1, R = I; are given in Tables. 11, 12 and 13. The split 
weighting case Is also considered and the corresponding closed loop 
poles are given in Table 14, The maximum actuator forces for the 
various cases described above are given in Tables 15, 16, and 17. 

To have an understanding of the relation between the actual coordi- 
nates and the modal coordinates (as the complete eigenvector matrix is 
not available) certain initial conditions are assumed for the modal 
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coordinates (q^. !■!, — 13) and the corresponding displacements are 
calculated. The nodes that experience maximum displacements together 
with the values of the displacements are shown In Table 18. It can be 
seen that nodes 100 and 101 have experienced the maximum displacements 
for the various Initial conditions assumed for the modal coordinates. 

The general conclusions for the various cases described so far are 
enumerated below. 

1. As the number of actuators Is Increased, for a given set of 
Initial conditions, the maximum force amplitude required for 
any actuator Is, In general, reduced. 

2. In general, by Increasing all elements In the Q matrix In the 
same relative amount, transient performance In the least damped 
modes Is Improved, but at the expense of a larger control effort. 

3. Split weighting of appropriate elements Ir. the Q matrix can 
result In Improved transient performance without unduly Increasing 
the control effort. 

4. System transient characteristics are degraded when the hoop 
mounted actuator Is not Included. Degradation is more pronounced 
as the number of (remaining) actuators Is reduced. 

5. For most cases simulated surface torsion and feed mast torsion 
are among modes with largest time constants. 

6. For the same Initial displacements in each of the 13 modal 
coordinates. It Is seen that the maximum translational dis- 
placements are realized at nodes 100, 101 (at the ends of the 
lower solar panels) . 


Table 7 
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COMPARISON OF CLOSEDLOOP POLES HOOP/COLUMM 13 ACTUATORS-13 MODES 

13 Actuators Poles 0=10001, R=I 13 Actuators Poles 0=10001 

Except 0(7,7)=Q( 10, 10)=:0( 12, 12) 
=0(20, 20) =0(22, 22) =0(25, 25)= 
100001, R=I 


( Real ) 

jw( Itnagi nary) 

(Real) 

Jw( Imaginary) 

-0.0030 

5.5505 

-0.0120 

5.5905 

-0.0080 

-5.5905 

-0.0120 

-5.5905 

-0.7045 

0.7788 

-0.4179 

0.4545 

-0.7045 

-0.7788 

-0.4179 

-0.4545 

-0.4179 

0.4544 

-1.0922 

1.5579 

-0.4179 

-0.4544 

-1.0922 

-1 .5579 

-2.1916 

1.6704 

-2.1916 

1 .6704 

-2.1916 

-1.6704 

-2. 191 6 

-1 .6704 

-0.9306 

2.6637 

-3.6157 

1.7991 

-0.9806 

-2.6637 

-3.6157 

-1.7991 

-3.6151 

1 .7994 

-1.0055 

0 

-3.6151 

-1.7994 

-1.0120 

0 

-1 .0058 

0 

-1.0250 

0 

-1 .0117 

0 

-1.0260 

0 

-1 .0247 

0 

-1.0864 

0 

-1 .0260 

0 

-1.0987 

0 

-1.0868 

0 

-1.1993 

0 

-1.0987 

0 

-1.290^ 

0 

-1 .2914 

0 

-1.5344 

0 

-1.5808 

0 

-2.5898 

0 

-3.9913 

0 

-3.5914 

0 

-5.4766 

0 

-5.4366 

0 

-6.2706 

0 

-6.3240 

0 

-10.0190 

0 

-10.06 17 

0 

-12.4240 

0 

-25.4183 

0 

-30.1860 

0 

-30.1864 

0 
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Table 8 
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COMPARISON OF MAXIMUM ACTUATOR FORCES HOOP/COLUMN 13 ACTUATORS- 

13 MODES 

Max. Actuator Forces 13 Actuators 13 Modes 13 Actuators 13 Modes 
(Amp.)-Pounds Qi ( 0) =0 . 0 1 , i = 1 , . . , 1 3 Qi ( 0 )= 0.01 , 1=1 , . . , 13 

0=10001, Rsl 0=10001, Except 

Q(7,7)=Q( 10,10)= 

Q( 12, 12)=0(20,20)= 

Q(22,22)=0(25,25)= 


f(1) 

0.3500 

100001, R=I 
1.7450 

f(2) 

0.0569 

0.0061 

f(3) 

0.3000 

1 .7330 

f(4) 

0.0569 

0.0061 

f(5) 

1.3027 

1.3029 

f(6) 

0.2833 

0.7333 

f(7) 

1 .2761 

2.3530 

f(8) (ft-lb) 

0.0123 

0.0317 

f(9) 

0.2369 

0.2277 

f( 10) 

0.1570 

0.1693 

f(11) 

0.M061 

0.4173 

f(12) 

0.3521 

0.3521 

f(13) 

0.0688 

0.0978 
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COMPARISON OF CL JSEDLOOP POLES HOOP/COLUMN 12/13 ACTUATORS 

Q=1000I,R=I 


13 Mode Model 13 Mode Model W.O Actuator? 

13 Actuator- Poles 12 Actuator-Poles 


(Real) 

jw( Imaginary) 

(Real) 

jw( Imaginary) 

-0.008 

5.5905 

-0.0024 

5.5905 

-0.003 

-5.5905 

-0.0024 

-5.590,. 

-0.7045 

0.7788 

-0.0162 

3.1813 

-0.7045 

-0.7788 

-0.0162 

-3.1813 

-0.4179 

0.4544 

-0.2654 

0.5653 

-0.4179 

-0.4544 

-0.2654 

-0.5658 

-2.1916 

1.6704 

-0.4179 

1 .6704 

-2.1916 

-1 ,6704 

• 0.4179 

-1.6704 

-0.9806 

2.6637 

-0.9173 

1.7994 

-0.9806 

-2.6637 

-0.9173 

-1.7994 

-3.6151 

1 .7994 

-2.1916 

0 

-3.6151 

-1.7994 

-2.1916 

0 

-1 .0053 

0 

-3.6153 

0 

-1.0117 

0 

-3.6153 

0 

-1 .0247 

0 

-1 .0117 

0 

-1 .0260 

0 

■ -1 .0247 

0 

-1.0863 

0 

-1 .0258 

0 

-1.0987 

0 

-1.0359 

0 

-1 .2914 

0 

-1 .0987 

0 

-1 .5808 

0 

-1.2785 

0 

-3.9913 

0 

-1 .6052 

0 

-5.47f6 

0 

-3.9913 

0 

-6,2705 

0 

-5.iJ867 

0 

-10.019 

0 

-6.3275 

0 

-12.424 

0 

-10.0673 

0 

-30.186 

0 

-30.1360 

0 


j 
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Table 10 
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COMPARISON OF ACTUATOR FORCES VARYING THE ACTUATORS WITH SEVENTH 
ACTUATOR REMOVED 


0=10001, R=I 


Max. Act. 
Forces(Amp) lbs 

7 Actuators 

f(1) 

1.4019 

f(2) 

0.4491 

f(3) 

1 .1369 

f(4) 

0.3411 

f(5) 

0.3612 

f(6) 

0.1209 

f(7) 

f(3) lb. ft 

f(9) 

f(10) 

f(11) 

f(12) 

f(13) 

0.2296 


10 Actuators 

12 Actuators 

1 .9990 

1.5187 

0.5261 

0.0176 

1.3024 

1.5293 

0.6535 

0.0175 

0.0636 

1.3024 

0.2369 

0.5597 

O'. 2047 

0.0423 

0.3603 

0.2279 

0.3520 

0.1695 

0.1072 

0.3636 

0.3520 

0.0739 


Table 11 
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COMPARISON OF CLOSEDLOOP POLES HOOP/COLUMN SYSTEM WITH SEVENTH 
ACTUATOR ON HOOP REMOVED- 13 MODES 0=10001, R=I 

7 Actuators 10 Actuators 12 Actuators 

(real) jw( Imaginary ) (Real) Jw( Imaginary) (Real) jw( Imaginary) 


0.0004 5.5905 

0.0004 -5.5905 
0.0007 3.1813 

0.0007 -3.1813 
0.0063 0.7467 

0.0063 -0.7467 
0.0434 0.0433 

0.0434 -0.0433 
0.4343 0.4053 

0.4343 -0.4053 
0.4690 1.3656 

0.4690 -1.3656 
0.5243 0.5476 

0.5243 -0.5476 
3.5693 1.367^ 

3.5693 -1.8674 
1.0180 0 
1 .0287 0 

1.0681 0 
1.2746 0 

1.6147 0 

2.4285 0 

3.4220 0 

6.0663 0 

10.656 0 
18.141 0 


-0.0016 

5.5905 

-0.0016 

-5.5905 

-0.0010 

3.1813 

-0.0010 

-3.1813 

-0.3084 

0.6657 

-0.3084 

-0.6657 

-0.4179 

0.4544 

-0.4179 

-0.4544 

-0.4505 

0.4712 

-0.4505 

-0.4712 

-2. 1<'16 

1 .6704 

-2.1916 

-1 .6704 

-3.6515 

1 .8369 

-3.6515 

-1.8369 

-1.0150 

0 

-1.0247 

0 

-1.0365 

0 

-1 .0907 

0 

-1 .0987 

0 

-1 .2839 

0 

-1 .5947 

0 

-3.9913 

0 

-3.9989 

0 

-6.1550 

0 

-10.661 

0 

-30.186 

0 


-0.0024 5.5905 

-0.0024 -5.5905 
-0.0162 3.1813 

-0.0162 -3.1813 
-0.2654 0.5673 

-0.2654 -0.5678 
-0.4179 0.4544 

-0.4179 -0.4544 

-0.9173 0.5953 

-0.9173 -0.5968 
-2.1916 1.6704 

-2.1916 -1.6704 
-3.6153 1.7994 

-3.6153 -1 .7994 
-1.0117 0 
-1.0247 0 

-1.0258 0 

-1.0858 0 

-1.0987 0 

-1.2785 0 

-1.6052 0 

-3.9913 0 

-5.4867 0 

-6.3275 0 

-10.067 0 

- 30.186 0 
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Table 12 


COMPARISON OF CLOSEDLOOP POLES HOOP/COLUMN SYSTEM WITH SEVENTH 
ACTUATOR REMOVED 13 MODES 0=100001, R=I 


7 

( Real) 

Actuators 
jw( Imaginary) 

( Real ) 

10 Actuator 5 
jw( Imaginary) 

12 , 
(Real) 

Actuators 
jw( Imagi 

-0.0005 

5.5905 

-0.0038 

5,5905 

-0.0062 

5.5905 

-0.0005 

-5.5905 

-0.0038 

-5.5905 

-0.0062 

-5.5905 

-0.0008 

3.1813 

-0.0230 

3.1812 

-0.0329 

3.1809 

-0.0008 

-3.1813 

-0.0230 

-3.1812 

-0.0329 

-3.1809 

-0.0113 

0.7467 

-0.2550 

0.5623 

-0.2497 

0.5597 

-0.0113 

-0.7467 

-0.2550 

-0.5623 

-0.2497 

-0.5597 

-0.0774 

0.7693 

-0.5012 

0.5193 

-0.5012 

0.5193 

-0.0774 

-0.7693 

-0.5012 

-0.5193 

-0.5012 

-0.5193 

-0.4855 

0.5466 

-1.3077 

0.4208 

-1.0011 

0 

-0.4855 

-0.5466 

-1.3077 

-0.4203 

-1.0023 

0 

-0.60S7 

0.1463 

-1 .0014 

0 

-1 .0023 

0 

-0.6087 

-0.1463 

-1 ,0023 

0 

-1 .0081 

0 

-0.6476 

0.6713 

-1 .0033 

0 

-1.0036 

0 

-0.6476 

-0.6713 

-1.0085 

0 

-1.0124 

0 

-1 .0017 

0 

-1.0036 

0 

-1.0549 

0 

-1.0026 

0 

-1 .0125 

0 

-1.3041 

0 

-1 .0062 

0 

-1.3041 

0 

-1.3257 

0 

-1.0124 

0 

-1.3215 

0 

-3.7535 

0 

-1.0924 

0 

-6.3906 

0 

-6.4064 

0 

-1.2619 

0 

-8.7613 

0 

-8.7613 

0 

-6.4271 

0 

-11.895 

0 

-12.655 

0 

-12.387 

0 

-13.193 

0 

-13.193 

0 

-13.099 

0 

-14.993 

0 

-19.119 

0 

-27.282 

0 

-27.292 

0 

-27.313 

0 

-34,811 

0 

-34.828 

0 

-34.853 

0 

-62.659 

0 

-98.551 

0 

-98.551 

0 
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COMPARISON OF CLOSEDLOOP POLES HOOP/COLUMN SYSTEM WITH SEVENTH 
ACTUATOR ON HOOP REMOVED 13 MODES Q=100000I, R=I 


7 

Actuators 

10 

Actuators 

12 Actuators 

( Real) 

jw( Imaginary) 

( Real) 

jw( Imaginary) 

(Real) 

jw( Imaginary) 

-0.0007 

5.5905 

-0.0079 

5.5905 

-0.0096 

5.5905 

-0.0007 

-5.5905 

-0.0079 

-5.5905 

-0.0096 

-5.5905 

-0.0011 

3.1813 

-0.0369 

3.1807 

-0.0399 

3.1806 

-0.0011 

-3.1813 

-0.0369 

-3.1307 

-0.0399 

-3.1306 

-0.0314 

0.7468 

-0.2488 

0.5592 

-0.2433 

0.5590 

-0.0314 

-0.7468 

-0.2488 

-0.5592 

-0.2483 

-0.5590 

-0.1383 

0.1358 

-0.5429 

0.5572 

-0.5249 

0.5572 

-0.1383 

-0.1358 

-0.5429 

-0.5572 

-0.5249 

- 0.5572 

-0.4809 

0.5699 

-1.0001 

0 

-1 .0001 

0 

-0.4809 

-0.5699 

-1 .0002 

0 

-1 .0002 

0 

-0.6212 

1.4837 

-1 .0003 

0 

-1 .0002 

0 

-0.6212 

-1.4837 

-1 .0008 

0 

-1.0008 

0 

-0.6585 

0.7347 

-1.0009 

0 

-1.0009 

0 

-0.6585 

-0.7347 

-1.0012 

0 

-1.0012 

0 

-1 .0000 

0 

-1 .0211 

0 

-1.0047 

0 

-1 .0002 

0 

-1 .0293 

0 

-1 .0293 

0 

-1 .0006 

0 

-1.0357 

0 

- 1.0349 

0 

-1.0013 

0 

-5.9089 

0 

- 12.513 

0 

-1 .0085 

0 

-20.433 

0 

-20.487 

0 

-1 .0232 

0 

-28.795 

0 

-28.795 

0 

-20.551 

0 

-37.884 

0 

-40.434 

0 

-39.367 

0 

-41.900 

0 

-41.900 

0 

-41.914 

0 

-47.751 

0 

-60.609 

0 

-88.147 

0 

-88.174 

0 

-88.248 

0 

-110.47 

0 

-110.53 

0 

-110.60 

0 

-199.66 

0 

-312.59 

0 

-312.59 

0 
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COMPARISON OF CLOSEDLOOP POLES HOOP/COLUMM SYSTEM WITH SEVENTH 
ACTUATOR REMOVED 13 MODES Q=1000I EXCEPT Q(7 , 7 ) =Q( 10 , 10) =Q( 1 2 , 1 2) 
=Q( 20,20) =0(22,22) =0(25,25) =100001 , R=I 


' 7 

Actuators 

10 

1 Actuators 

12 Actuators 

1 (Real) 

jw( Imaginary) 

(Real) 

jwCIroaginary) 

(Real) 

jw( Imaginary) 

{ - 0.0013 

5.5905 

-0.0050 

5.5905 

-0.0068 

5.5905 

j -0.0013 

-5.5905 

- 0.0050 

-5.5905 

-0.0063 

-5.5905 

1 -0.0024 

3.1813 

-0.0277 

3.1811 

-0.0381 

3.1303 

-0.0024 

-3.1813 

-0.0277 

-3.181 1 

-0.0381 

-3.1808 

-0.0200 

0.7467 

-0.2389 

0.3015 

-0.2367 

0.3100 

-0.0200 

-0.7467 

-0.2389 

- 0.3015 

-0.2367 

- 0,3100 

-0.0434 

0.0433 

-0.4179 

0.4544 

-0.4179 

0.4544 

-0.0434 

-0.0433 

-0.4179 

-0.4544 

-0.4179 

- 0.4544 

-0.4343 

0.4053 

-0.9383 

0.6871 

-2.1916 

1 .6704 

- 0.4343 

-0.4053 

-0.9383 

-0.6871 

-2.1916 

-1.6704 

-0.4690 

1 .3656 

-2.1916 

1.6704 

-3.6158 

1 .7991 

-0.4690 

-1.3656 

- 2.1916 

-1.6704 

-3.6158 

-1 .7991 

-0.5240 

0.5478 

-3.6401 

1.8295 

-1.0017 

0 

-0.5240 

-0.5478 

-3.6401 

-1 .8295 

-1.0247 

0 

i -3.5693 

1 .8675 

-1 .0150 

0 

-1 .0258 

0 

-3.5693 

-1 .8675 

-1 .0248 

0 

-1 .0853 

0 

-1 .0180 

0 

-1.0349 

0 

-1.0987 

0 

-1 .0287 

0 

-1 .0879 

0 

-1 .1634 

0 

-1.0680 

0 

-1.0987 

0 

-1 .2852 

0 

- 1.2750 

0 

-1 .2831 

0 

-1.5958 

0 

-1.6150 

0 

-1.5966 

0 

-2.6473 

0 

-2.4280 

0 

-3.9914 

0 

'3.9913 

0 

-3.4220 

0 

- 4.3071 

0 

-5.4839 

0 

1 -6.0669 

0 

-6.1674 

0 

-6.3652 

0 

\ -10.066 

0 

-10.661 

0 

-10.673 

0 

1 -18.141 

r 

0 

-30.186 

0 

-30.186 

0 


I 
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Table 15 


COMPARISON OF ACTUATOR FORCES VARYING THE ACTUATORS WITH SEVENTH 
ACTUATOR REMOVED 


Q=10000I, R=I 


Max, Act. 7 Actuators 

Forces(Amp) lbs 

f(l) 4.0402 

f(2) 3.3862 

f(3) 4.3702 

f<4> 2.3451 

f(5) 1.1247 

f(6) 0.2136 

f(7) 

f(8) lb. ft 1.2184 

f(9) 

f(10) 

f(11) 

f(12) 

f(13) 


10 Actuators 

12 Actuators 

4.4356 

3.5095 

3.5324 

0.0958 

3.6011 

3.3440 

1 .4979 

0.0953 

0.2539 

3.6010 

0.3500 

1.2351 

0.b595 

u.l'721 

1.1955 

0.3503 

0.9433 

0,6543 

0.3156 

1 . 1675 
0.9437 
0.2844 




Table 16 
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COMPARISON OF ACTUATOR FORCES VARYING THE ACTUATORS WITH SEVENTH 
ACTUATOR REMOVED 


0=1000001, R=I 


Max. Act. 
Forces(Amp) lbs 

7 Actuators 

10 Actuators 

12 Actuators 

f(1) 

9.6904 

3.4849 

6.5931 

f(2) 

20.7948 

8.6722 

4.5791 

f(3) 

14.6512 

1 .1022 

5.6322 

f(4) 

8.8780 

3.5757 

0.4579 

f(5) 

3.5327 

1 .4109 

1.1023 

f(6) 

1 .2470 

3.3063 

3.1263 

f(7) 

f(8) lb. ft 

4. 1b67 

2.7808 

1.2149 

f(9) 


3.9409 

3.3061 

f(10) 


2.4642 

2.8507 

f(11) 


1 .0831 

3.9445 

f(12) 



2.4643 

f(13) 

Table 

17 

1 .0771 

0=10001, EXCEPT 
Q(25,25)=10000I, 

Q(7,7)-Q(10,10) 
, R = I 

=Q(12,12)=Q(20, 

,20)=Q(22,22)= 

Max. Act. 
Forces(Amp) lbs 

7 Actuators 

10 Actuators 

12 Actuators 

f(1) 

3.0369 

4.9372 

3.5906 

f(2) 

1 .2538 

1.1888 

0.0696 

f(3) 

1.1868 

1 .3020 

3.6637 

f(4) 

1.4440 

1 .5255 

0.0696 

f(5) 

0.3660 

0.2012 

1.3022 

f(6) 

0.1209 

0.2367 

1.2361 

f(7) 




f(8) lb. ft 

0.5030 

0.2294 

0.1461 

f(9) 


0.3169 

0.2369 

f(10) 


0.3520 

0.1751 

f(11) 


0.2883 

0.3452 

f(12) 



0.3520 

f(13) 



0.1141 
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RELATIONSHIP BETWEEN DISPLACEMENT IN THE MODAL COORDINATES AND 
MAXIMUM TRANSLATIONAL DISPLACEMENT (AT NODE 101) 


Node 


Inti. Condition 

Max. translatio- 



Qi(0) 

nal Disp 

.-Inches 

(100) 

101 

0.01 

(0.2840) 

0.3009 

(100) 

101 

0.02 

(0.5630) 

0.6019 

(100) 

101 

0.03 

(0.8519) 

0.9029 

(100) 

101 

0.04 

(1.1359) 

1 .2039 

(100) 

101 

0.05 

(1.4199) 

1 .5048 

(100) 

101 

0.06 

(1.7039) 

1.8058 

(100) 

101 

0.07 

(1.9878) 

2.1068 

(100) 

101 

0.08 

(2.2718) 

2.4078 

(100) 

101 

0.09 

(2.5553) 

2.7037 

(100) 

101 

0.10 

(2.8398) 

3.0097 
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m. ORIENTATION AND SHAPE-CONTROL OF AN 

ORBITING FLEXIBLE BEAM UNDER THE INFLUENCE 
OF SOLAR RADIATION PRESSURE 


In this paper, the uncontrolled and controlled dyna- 
mics of a thin flexible beam in orbit and in the 
presence of, solar radiation disturbance are analyzed. 
A beam nominally oriented along (1) the local hori- 
zontal and carrying a gimballed rigid dumbbell for 
gravity stabilization, and (i1) a beam nominally 
oriented along the local vertical are considered. 

The uncontrolled dynamics of the beam in the pres- 
ence of the solar radiation pressure disturbance 
shows the excitation of the rigid pitch mode. The 
control laws previously designed for the case where 
the environmental effects were neglected, are found 
to be inadequate to control the shape and orienta- 
tion of very flexible beams that are exposed to 
solar radiation disturbances. The control laws 
and the gain parameters are reevaluated for both 
cases of nominal, beam orientations; this results, in 
general, in Increased robustness, of the closed-loop 
system. Methods of obtaining a robust control sys- 
tem in the presence of environmental perturbations 
are discussed. 


INTRODUCTION 

At the operational altitudes of proposed future missions involving 
large space structures, the principal environmental disturbance will, 
in general, be due to solar radiation pressure. Solar pressure torques, 
induced on the structure, are dependent on the surface properties and 
the surface geometry.! Expressions for the induced solar pressure 


torques on a free-free flexible beam were developed earlier.^ The un- 
controlled dynamics of a thin flexible beam in orbit under the in- 
fluence of solar radiation pressure forces were considered for two cases 
of beam orientations, namely, (1) nominally oriented along the local 
vertical and (11) nominally oriented along the local horizontal and gra- 
vity-stabilized using a rigid dumbbell attached at the center of mass 
of the beam. The uncontrolled and controlled dynamics for the two cases 
of nominal beam orientations have also been considered previously but 
in the absence of environmental disturbances.^*^ 

The objective of the present study Is to evaluate the validity of the 
previously developed control laws, which were obtained disregarding 
the environmental effects, for controlling the shape and orientation of 
a beam actually under the Influence of solar radiation disturbance. It 
Is proposed to evaluate the robustness of the previously developed con- 
trol laws for the beam under the Influence of the solar pressure dis- 
turbances by considering a parametric study of the controlled dynamics 
for various Initial conditions. Particular attention will be given to 
highly flexible structures under larger Initial displacements. The 
feedback gains will be modified, where required, to achieve satisfactory 
transient performance while minimizing control force effort. Where 
modification of gain values are found not effective, modification of 
control laws and reselection of actuators and their locations will be 
considered. 

THE MATHEMATICAL MODELLING OF A FLEXIBLE ORBITING BEAM 

Dumbbell Stabilized Flexible Beam in Orbit Nominally Oriented Along 
the Local Horizontal 


A flexible orbiting beam nominally oriented along the local horizontal 
represents a gravitationally unstable configuration. A passive sta- 
bilization of the beam can be obtained by using a rigid dumbbell with 
proper moment of inertia. In Ref. 3, the equations of motion for a 
beam with a dumbbell assumed to be attached at the center of mass of 
the beam (Fig. la) through a spring loaded hinge and with viscous rota- 
tional damping have been developed. It Is assumed that the dumbbell 
mass is concentrated at the tips, that the viscous force at the hinge 
Is linear, and that all displacements and deformations occur within the 
orbital plane. With the usual assumptions of small pitch amplitude and 
dumbbell oscillations and flexural deformations, the linearized equa- 
tions of motion are obtained as^. 


e«+c0'+{F-3)0-ca'-Fa + E(ce;;+ke^)C^"^=» Cy/OyW^ (1) 

- a"+c^to'+(c^F+3)o-CiF0'-c.jF0- £(Fe^+Fejj)c-jC^”^ “ ® 

= E /M (3) 

n' n c 
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where “ 1,2,...) and = n'^h modal mass 

z = undeformed beam length 

6 = pitch angle 

» orbital angular velocity 
Cy = external torque about the pitch axis 

Jy » beam pitch axis moment of inertia 

= nth flexural modal amplitude 
S “ '^o'i‘l'*'"6'^sionalized nth flexural modal amplitude 

V » k/JyW^; c = c/JyW^ 

k = torsional restoring spring constant at the hinge 
c = viscous damping coefficient 

a = angle between the dumbbell axis and the local vertical 

r(n) . 

^ ^ X.0 ' 

s beam shape function of the nfh transverse mode 
c^ = “ pitch moment of inertia of the dumbbell 

= w^/wj. j Wyj = n^*^ natural frequency 
( )' = d/dx where t“ o)^t 
t * time 

The following observations can be made from a study of Eqs. (1),(2) 
and (3) in the absence of external forces on the system: (aj the 

pitch motion of the beam, the dumbbell motion (a), and the elastic 
motion of the beam (e ) are all coupled to each other; (b) within the 
linear range the elastic modes for which 0 (the symmetric modes), 
are completely independent of the pitch and dumbbell motions. Further- 
more these modes do not influence either the pitch or the dumbbell mo- 
tion. 

A Flexible Beam in Orbit Nominally Oriented Along the Local \fertical 

The equations of motion for a thin uniform flexible beam in orbit with 
its axis nominally along the local horizontal (Fig. lb) are developed 
in Ref, 3 and represent a gravitationally stable system. With the 
assumptions of small transverse in-plane vibrations and small ampli- 
tude pitch oscillations of the beam, the linearized equations of 
motion are obtained as^» 

9”+30 . original PAGE 13 W 

of poor quality 
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Active control is required to maintain the pitch orientation of the 
beam and to damp out the modal oscillations in the presence of the 
environmental disturbances. A point actuator located at one end of 
the beam is considered for shape and orientation control of both the 
dumbbell stabilized beam and the beam oriented along the local verti- 
cal. The development of the forcing terms on the right hand side 
of Eqs. (l)-(5)» due to the solar radiation pressure is reviewed in 
the next section. 

SOLAR RADIATION DISTURBANCE MODEL 

A detailed study of the effect of solar radiation disturbance on the 
uncontrolled dynamics of a flexible beam appears in Ref. 2 and a 
summary of the disturbance model obtained is given below. Fig. 2 
shows the geometry Of reflection of a flexible beam vibrating in one 
of the free beam modes and exposed^to solar radiation pressure. The 
incidence solar radiation vector, t, and the unit normal to the sur- 
face, n, are assumed^to be__in the plane of the transverse vibration 
of the beam. Then, t and n" can be expressed as 

T “ a^jT + Cq1< (6) 

n = (((»'T-F)//Tfi^ (7) 

where, <|)' * d(j>/dC, <}> is the free-free beam shape function and 5 is 
the nondimensionalized longitudinal coordinate of the beam. 

The Solar radiation force, fg and moment, Ng, acting on a completely 
absorbing surface is given by^ 

* h^T / T*n ds (8) 

5 

N. =“h 7 X / R(T*rT) ds (9) 

a 0 s 

wiiere h = 4.64x10"® Nt/m^ is a constant for earth orbiting space- 
craft. °The integration over the area, s, is bounded by the condition 

( 10 ) 

sur- 


( 11 ) 

( 12 ) 

of 


T*n ^ 0 

The corresponding force and moments for a completely reflecting 
face are given uy 

n 1 * 2h_ / rT X R(T*rr)^ ds 
I ® s 

where IT is the position vector of ds with respect to the center 
mass of the beam. For a surface with an arbitrary reflection 
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(13) 

(14) 


coefficient, c^, the force and moment expressions become': 

%r * fa *''r(%-''^a) 

\r ■ Na^MV^a) 


Only the moment expressions are of importance in analyzing the atti- 
tude motion of the beam. The expressions for the moments per unit 
width of the beam are developed -using Eqs. (6) and (7) In Eqs. (9) 
and (12) as2: 



"*^o®o*'o^o symmetric modes) 

0 (for asyronetrlc modes) 



,1 (ao’f’-'o'^ 


{ <(,'4, + (5- 





(15) 

(16) 


The numerical Integration of Eq. (16) has shown that the moment due 
to any asymmetric mode of the beam Is extremely small as compared 
with the magnitudes for the symmetric modes and the moments, 
due to the symmetric modes can be expressed 1n the form (similar to 
Eq. (15)), 


Nr “ ®Q^o^o synmetrlc modes) (17) 

Eq. (17) Is valid for all symmetric modes and for small deflections 
of the beam. Eqs. (15) and (17) are used In Eq. (14) to obtain the 
total moment due to solar radiation pressure on a beam with a coef- 
ficient of reflectivity, c-. Further, the deflection at one end of 
the vibrating beam for a given mode can be written as 

«n(t) =» e„(t) l 

Thus, the total moment on the beam due to any one of the symmetric 
modes of the beam can be expressed as 


N(t) . c„(t) 1. a^jC^ (18) 

where a. = sin 6, and c„ = cos 6, and 3, is the solar Incidence angle 
given by° ’ ® ’ 


0^(t) = Wj.t + e(t) + 0 ^( 0 ) (19) 

Where ® orbital angular velocity 
0(t) * pitch angle of the beam. 

ORIGINAL PAGE 13 
OF POOR QUALmr 


3.5 


The effect of the disturbance on the generic mode is obtained b^ 
evaluating the integral® 




r 


T ds 


( 20 ) 


where, F is the external force due to solar pressure. Eqs. (8) and 
(11) are substituted into Eq. (20) and the resulting integrals are « 
evaluated under small amplitude approximations with the result that. 


E » E 
na n. 


and 


'a 

» E 


ttg r Oj. n^ 


( 21 ) 


Thus, the principal effect of solar radiation pressure on the flex- 
ible beam is to produce a net moment about the pitch axis of the beam.. 
For a beam of length 100m and c« = 0.5, the maximum pitch moment per 
unit deflection, Nm, in Eq. (18) is obtained using Eqs. (14), (15) 
and (16) as,^ 

» 2.23x10"^ + 0.5(9.4xl0"^-2. 23x10"^) 

* 1.58x10'^ Nt-m 



calculated here is for a beam with unit width). Eq. (18) is 
ready to be used as the pitch forcing function in Eqs. (1) and (4). 


CONTROL OF THE ORBITING BEAM SYSTEM 


For applications of control theory the sets of second order equations, 
Eqs. (l)r(3), and (4) and (5) are transformed into a state vector 
form given by 

X = AX + BU + D (22) 

where X is the state vector of order (2m), m is the total number of 
rigid and flexible modes in the system. The plant matrix has dimen- 
sionality (2mx2m), B, (2mxr), where r denotes the number of actuators 
and D is a vector of order 2m denoting the disturbance and is a 
function of both time and the state of the system. For the dumbbell 
stabilized beam with two flexible modes included in the model, the 
order of the system will be 8, and for the beam oriented along the 
local vertical the order of the system v/ith the two flexible r,. »des, 
will be 6. 


A single point actuator located at one end of the beam was considered 
for both cases of the beam's nominal orientation in Ref. 4. An ideal 
feedback law of the form U = -KX is assumed to obtain a control law 
and the corresponding gain values, K. In the following sections the 
control problem is considered separately for the two cases of the 
Learn orientation. 


Control of the Dumbb^l Stabl 1 1 zed FI ex i ble Beam 


The disturbance due to the solar radiation effects on the dumbbell Is 
neglected in the following study, since the dumbbell Is assumed to 
be rigid. The solar radiation disturbance vector, D, is then given 
by 


D * CO, 0.'' n, 0,0,0]^ 

(23) 

/ 

N » e.j£ Njjj sin 6,j cos 6^ (from Eq. 18) 

(24) 


In the abs«jce of both the control and the disturbance, an Initial 
displacement In the first modal amplitude would correspond to simple 
harmonic motion In the first mode only. With the Inclusion of the 
solar radiation disturbance (given by Eq. (23)) into the model, it is 
seen that the pitch motion is also excited. As an example. Fig. 3 
illustrates the response of the beam to an Initial displacement of 
£,(0) = O.O"’ in thp presence of the solar radiation disturbance where 
the fundamer.t >1 flt.- jral frequency is assumed to be ten times the 
orbital rate. The pitch motion also Induces the dumbbell motion and 
oscillations in the second mode through coupling. The Induced pitch 
amplitude is nearly 2®. 

The method adopted hereto obtain control laws and the gain parameters 
for which an acceptable transient response would result In the pre- 
sence of the disturbance is to first consider the control problem by 
Ignoring the disturbance, D, in Eq. (22,). The transient response of the 
system with the disturbance is obtained next, and then the gain para- 
meters and, if necessary, also the form of the control laws will be 
modified to obain an acceptable system response under the action of 
the disturbances. 

A set of gain parameters was obtained earlier by examining the charac- 
teristic equation of the dumbbell stabilized beam with one actuator 
at one end of the beam, for which the disturbances were completely 
Ignored. The control law obtained was of the form 

F^ = -0.070 -O.OOSe^ -0.03c2 (25) 

The transient responses with this control law for the dumbbell sta- 
bilized beam with and without the disturbances Included in the system 
are shown in Fig. 4 for an initial displacement of O.OU in the first 
mode only. The fundamental frequency of the beam Is assumed to be 
ten times the orbital frequency. The effect of the solar radiation 
disturbance is to increase the induced pitch amplitude slightly as 
shown in Fig. 4. The peak control force required Is also increased 
by a small amount, but still. is of the order of 10“^ Nt. only. 

the control law for this case is seen to be acceptable under 
the expected solar radiation disturbance on the beam.: 


To study the robustness of this control system, an extreme case of 
the solar radiation disturbance is considered with an initial dis- 
placement of O.U In the first mode for a beam whose fundamental fre- 
quency is 7 times the orbital frequency. For this case of the beam 
and based on the previous control law, Eq. (25), the transient re- 
sponse of the system Is shown In Fig. 5 with and without the solar 
radiation disturbance. It Is clearly seen that the dumbbell sta- 
bilized beam whose pitch motion Is well controlled in the absence 
of the solar radiation disturbance, experiences extremely large ampli- 
tudes of pitch oscillations (+150) due to the solar radiation dis- 
turbance. The transient performance of the system in all modes Is 
Improved simply by Increasing the gain value corresponding to the 
modal rate by ten times and also the gain proportional to the pitch 
angle from -0.07 to -0.3, (Fig. 6), for the same initial conditions 
and beam parameters as In Fig. 5. The large amplitudes In the pitch 
motion disappear as the first modal amplitude is quickly damped re- 
ducing the magnitude of the solar radiation disturbance. The peak 
control force now Increases by two orders of magnitude as compared 
with the cases in Fig. 3 and Fig. 4, but the peak force value of 
2xl0"2Nt. Is still a small value. Thus, the need to consider the 
solar radiation disturbance In designing a control system for a 
highly flexible system Is demonstrated. 

Control of a Flexible Beam Nominally Oriented Along the Local 
Vertical 


The same procedure used for the dumbbell stabilized beam is repeated 
here to obtain a robust control law and the gain parameters for use 
with the solar disturbance input. The control law, with the gain 
parameters selected to provide damping in pitch and tl.a first mode 
is given by 

Fg = -0.01 0 -0.01 (26) 

The second mode is not Included in this model, because the second 
mode Is decoupled from pitch and the first mode, even in the pre- 
sence of the solar radiation disturbance. ^ The system transient 
response with the solar radiation disturbance is shown in Fig. 7 
for e-](0) = O.OH and the fundamental frequency of the beam assumed 
to be ten times the orbital frequency, and with and without the 
application of control. In the absence of control, the magnitude of 
the Induced pitch oscillations (*0.25°) is small as compared with 
the case of the dumbbell stabilized beam (Fig. 3), because of the 
stabilizing gravity forces. With the application of the control, 
damping in both pitch and the modal amplitude results 

The same control law, (Eq. 25), applied to a more flexible beam 
(wi = 3.0) and a larger initial amplitude in the first mode 
(ei(0) = O.U) shows (Fig. 8) that the pitch amplitude overshoot 
is nearly 30°. An attempt was made to increase the feedback gain 
proportional to the first flexible mode, e-], as before, for the case 
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of the dumbbell stabilized beam. The transient responses (not shown) 
Indicated large overshoots In the pitch amplitude and a further In- 
crease in the feedback gain value made the system unstable. Similar 
attempts of varying the gain values corresponding to the pitch rate 
and providing the feedback proportional to the pitch angle failed 
to provide an acceptable transient response of the system. Two 
actuators were then assumed, one at the center of the beam (pro- 
viding control over the flexible mode only) and the other at one of 
the nodal point ( 4 » 0.776) of -the flexible mode to control the 
pitch motion only.*' The pitch and the flexible mode were then de- 
coupled except for the Influence of the flexible mode on the pitch 
motion through the solar radiation disturbance. The following control 
laws were selected for the two actuators to provide critical damping 
in both pitch and the flexible mode. 

fl = -0.1136 » -0.3e^ (27) 

The transient response for the same Initial conditions In Fig. 8 and 
under the new control law [Eq. (27)] is shown In Fig. 9. The pitch 
motion and the flexible modal oscillations are quickly damped out. 

The peak control force required Is of the order of 10"2Nt. for both 
the actuators - the s^e order of magnitude required In Fig. 8. 

CONCLUSIONS 

The uncontrolled and controlled dynamics of a thin flexible beam In 
orbit and in the presence of solar radiation disturbances are 
analyzed. Control laws and gain parameters are obtained to control: 

(I) the beam nominally oriented along the local horizontal and 
carrying a gimballed rigid dumbbell for gravity stabilization, and 

(II) the beam nominally oriented along the local vertical. The 
control laws previously developed Ignoring the environmental effects 
are found to be Inadequate to control the shape and orientation of 
very flexible beams that are exposed to solar radiation disturbances. 
In order to obtain a robust control system In the presence of envlro- 
mental perturbations, it is sometimes more desirable to Increase 

the number of actuators rather than simply modifying selected gain 
values. For further research related to the problem of minimizing 
the overshoot under large disturbances (in a qualitative' manner) the 
use of pole placement techniques and application of linear quadratic 
Gaussian method are suggested as a follow on effort, especially for 
large order systems. 
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Fig. 1(a) Dumbbell Stabilized Flexible Beam Nominally Oriented 
along the Local Horizontal. 



i 


Fig. 1(b) Oniform Flexible Beam Nominally Oriented along the 
Local Vertical. 
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Deflected Beam' 


I^deflected Beam 


Pig. 


. Geometry of Saflectioa of a Flexible Beam 
Exposed to Solar Badiatlon. 
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^lg« 3. ' Time Response of Dumbbell Stabilised Flexible Beam la the 
Presence of Solar Radiation Pressure, ■ 10.0), 


\ 
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Time Kesponse of Dijinbbell Stabilized Flexible Beam In Orbit 
with Modified Active Control and In the Presence of Solar 
Radiation Disturbance. 







without Control — — — With Control 





I.C.'s 9(0) - 0 Ej^(O) - 0.01 

F - - 0.01 6 - 0.01 e - 0 ), “ 10.0 

c 11 

Fig. 7 . Time Response of the Flexible Beam Nominally Oriented 
Along the I.ocal Vertical and in the Presence of Solar 
Radiation Pressure. 
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Time (orbits) 


e ^( 0 ) - 0.1 0 ( 0 ) - 0 - 3,0 

Time Response of Flexible Beam Nominally Oriented Along the 
Local Vertical and In Presence of Solar Radiation Pressure and 
With and Without Active Control. 
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IV. EFFECT OF SOLAR RADIATION PPilSSURE ON THE DYNAMICS OF A THIN 
HOMOGENEOUS SQUARE PLATE IN ORBIT 

IV.l Introduction 

Proposed future applications of large space structures require .ontrol 

of the shape and orientation of the structure In orbit. It has been shown 

previously^, considering a long, thin and uniform beam, that the principal 

environmental disturbance acting on these structures could be due to the 

solar radiation pressure. In the present work the dynamics of a more 

important basic structure, namely, a thin, homogeneous and flexible square 

plate exposed to solar radiation disturbance will be considered. The force 

2 

and moment expressions as given by Karymov will be used to obtain the 
expressions for solar radiation disturbing forces and moments acting on 
the free-free square plate in orbit. The dynamics of such a plate nominally 
oriented along the local vertical was considered earlier disregarding the 

3 

environmental disturbances. In the present study it is proposed to recon- 
sider the dynamics of the square plate nominally oriented along the local 
vertical with the solar radiation force and moment expressions included in 
the dynamic model. 

The mode shapes and the frequencies of the plate are obtained using 

4 

the finite element program, STRUDL. To obtain expressions for solar 
radiation forces and moments, it is convenient to express the mode shapes 
of the plate as a combination of the mode shapes of a free-free beam.^ 

The first five modes of the plate will be considered for study here. 

The plate is assumed to have only small transverse vibrations, so that the 

shadowing of the plate due to any deflected part of the plate can be 
neglected. The small deflection assumption also allows the superposition 
of the beam mode shapes in representing the mode shapes of the plate. 
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IV. 2 Solar Radiation Forces and Moments Acting on a Thin Homogeneous 
Flexible Square Plate 


Fig. 4.1 shows a square plate exposed to solar radiation. Let n 
denote the outward unit vector normal to the surface, ds. and let T 
be the unit vector in the direction of solar radiation denoted as 


T-al + bj+c k 
o o o 


(4.1) 


The direction cosines of T, namely, b^ and c^, can be expressed in 
termi, of the solar incidence angles,©^ and (defined in Fig. 4.1) 

dS 

a = sin 9. cos Ui. 

o i ^1 

b “ sin 9. sin ijJ. (4.2) 

O 1 

c “ cos 9. 

o i 

Then, the solar radiation force, F , and the moment, N , on a completely 

3 a 


absorbing surface are given by 


and 


F ■* h T / T*n ds 
o s _ 

N —h T X / R (T*n) ds 


(4.3) 

(4.4) 


o 

2 

where, h =• 4.64x10 Nt/m is a constant for near earth space structures. 


The Integration over the area, s, is bounded by the condition 


T*n ^ 0 (4,5) 

The force, F^, and moment, acting on a completely reflecting surface 

2 

can be developed as , 

F = -2h / n (T*n)^ ds (4.6) 

2 

and " 2h / nxr (x*n) ds (4.7) 

_ ® 

where, R is the position vector of ds with respect to the center of mass of 
the plate. 


4.2 


, ORIGINAL PAGE 19 

OF POOR QUALITY 

For a surface with an arbitrary reflection coefficient, Er, the force and 

2 

moment expressions become : 


F - F + e (F “F ) 
er a r' r a 

i - N + e (ii “N ) 

er a r r a*^ 


(4.8) 


The shape function of a rectangular plate can be represented as a 
product of the two beam functions given by^ (considering only the trans- 
verse vibration) , 

6 and are the free-free beam shape functions given by 

0 (x) = a (slnQ X + sinhfl x) + (cosJi x + coshfi x) 
n n n n n n 


for n = 2,3, 4 . . . . (4.10) 

where, O “ (cos^2 -coshfi )/(slnhfl -slnfi ) 
n n n n n 

and 9 (x) “ l-2x for n=l 
n 

“ a constant for n=0 

and 0„(^) = t (y) 

n n 

For a square plate, certain special modes which are combinations of the 
modes of a rectangular plate are shown to be existing.^ The frequency 
expressions for such modes are also given In Ref. 5. The first five 
modes of a square plate in which the second and third modes represent 
special combinations of "beam modes" (Fig. 4.2) are considered in the 
present study. 
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A unit normal to the syrface, n, la given by, 
n ■ a,l + b-j -I' c.k 
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*1* ■ “i*^ '•r 

- il;(de/dg)l + 9(d^/dn)T - k 

/ (\l> de/dC)2 + 6 (dijz/dn)^ + 1 
((CiH* C) are non-dimensional coordinates In the x,y,z directions, respectively.) 
The position vector, R, is represented as, 

R 


(C- |) 1 + (n - -J) J + 2 ^ 


(4.12) 

Eqs. (4.1), (4.11) and (4.12) are substituted into Eq. (4.4) and then the 
resulting integrals are evaluated to obtain the expression for the moment 
acting on a plate having a completely absorbing surface as, 

^a “ ^ ^%(Sj^-s^/2)-aQS3} J 

+ {a^(s2“S^/2) - b^ (sj^-s^/2)} k] 

where , 


(4.13) 


®1 " ^ 5s^dCdri 

s 

» / Cs d?dn 
^3 3 *^ 


f ns^d^dn 


s^ = / s^d^dn 


(a 'IL e_ 


di//, 


m 


- ‘^o^ 


o'^m d5 ’ “o ''n dn 
Aue Integrals to s^ can be evaluated analytically. The moment express- 
ions are obtained for the first five plate modes (Fig. 4.2) by evaluating 
^1 ®4 combinations of corresponding (m,n) modes and are given as, 

N 


. 2 


“ [b 1 + a j!I z, (for mode II) 

o 0 o o z 

« h^A^c [b 1 - a^jl z- (for mode III) 
o o o o J 

“ o, (for modes IV and V) 


(4,14) 
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where, ^2 ^3 deflections at one comer of the plate asso- 

ciated with the I, II and III modes, respectively. 

The moment due to solar radiation pressure, Nj., acting on a completely 
reflecting surface is obtained by substituting Eqs. (4.11), (4.12) and 
(4.1) in Eq. (4.7). The resulting integral is simplified to obtain the 
expression. 

Nr = 2h^ / [(a^C-a^ri )i-+ (a^? -aj^?)J +(aj^H *- 82 ? )k] s^d^dr) (4.15) 

s 

where, ? ' « ? - 0.5 and n' = H - 0.5 

Eq. (4.15) involvis complicated integrals and to find an analytical solu- 
tion is very difficult. Instead, a numerical evaluation of the integrals 
involving different modes are carried out and the results are shown in 
Fig. 4.3. The plate dimension is considered to be lOOmxlOOm and the 
deflection at the corner of the plate for each mode is assumed to be 
“ Z 2 = = 1.0m. Similar results are also shown in Fig. 4.3 for a 

plate having a completely absorbing surface and are obtained using Eq. 
(4.14). The solar incidence angle, 9^, is varied from 0 to 90°, with 

0, Only the first three modes give rise to appreciable moments for 
both completely absorbing and completely reflecting surfaces. The 

magnitudes of the moments are seen to be an order of magnitude higher 
-2 

(2x10 Nt-m) for a completely absorbing surface as compared with the case 

—3 

of a completely reflecting surface (10 Nt-m). The moments due to 
modes II and III, and for both completely reflecting and completely ab- 
sorbing surfac s, can be visualized as extensions of the result obtained 
for the case of the beam.^ For = 0° (incidence only in the z,x 
plane) the moment 01 a beam is shown to be aoout the y axis for any given 
mode of the beam. 
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But, for the same case, Interaction of the solar pressure on the plate 
vibrating in only the first mode, results in a moment about the x-axis. 
The reason for the latter phenomena can be explained by considering 
the force distribution for a completely absorbing surface (Fig. 4.4). 
There is a finite moment about the x-axls along the cross section 


y '-y * as shown in Fig. 4.4 because of the uniform variation in the 
angle of Incidence along y '-y All the cross sections must be 
similarly considered from 5= 0 to 1, and the net moment result will, 
in general, be non zero. A similar explanation can be given for the case 
of a completely reflecting surface. 

Based on the numerical results shown in Fig. 4.3, in which is 

varied from 0 to 90° (not shown) the moment expressions for a completely 

reflecting plate can be written as, 

h, c (a 1-b j) z, (for mode I) 
i o o o 1 




“ ^2*^0 ^2 mode II) 

" h„c (b i-a j) z^ (for mode III) 

A 0 O O 

where, " ■ 3.25x10“^ and h 2 = 1.09x10 

Eq. (4.16) is found to be valid for magnitudes of z^ to up to O.OIJI, 

The moments about the x,y and z axes are obtained by collecting the coeffi- 
cients of 1, j and k, respectively, from Eqs. (4.14) and (4.16) as, 


N 


ax 


N 


h3 {(a^/3)z^ + b^(z2+Z3)> 
-h^ {(b^/3)z^ ■^^o^^3"’*2^ ^ 


ay 

"a. " 


(4.17) 


where h- = h Jl c 
3 o o 
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( 4 , 18 ) 


IV. 3 Modal Forces Due to Solar Radiation Pressure 


The effect of the disturbance on the generic mode is obtained by 
evaluating the integral^ 

E ■ / z (x,y) k*dF (4.19) 

n 

where* F represents the force due to solar radiation pressure. Also, 

E - E + £ (E -E ) (4,20) 

n na r ' nr na 

Eq. (4.3) is substituted into Eq. (4.19) and after evaluating the re- 
sulting Integrals, E^^^ is found to be equal to zero for all modes of the 
plate. Eq. (4.6) is used in Eq. (4.19) to get 

E^^ - 2h^ / [z(a^aj^+b^bj^+c^)/(a^+b^+l)3 dCdn (4.21) 

The slopes, d0^/d^ and are assumed to be very small so that 

2 2 2 

a^^+bi+Ci K 1. Thus, the Integral in Eq. (4.21) can be easily evaluated 
to show that E^^ is also equf. to zero for all modes of the plate. Hence, 
the solar radiation pressure, does not give rise to any generic force. The 
results obtained can now be used in the d3mamic model of a flexible plate 
in the orbit. 
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IV. 4 Effect of Solar Radiation Pressure on a Plate Nominally Oriented 
Along the Local Vertical 

The major surface of the plate is assumed to be perpendicular to 
the orbital plane (Fig. 4.5). From the general formulation of Refs. 3 and 
5, the equations of motion of the structure are obtained, under the assump- 
tion that the transverse deformations are small compared to the charac- 
teristic length of the plate. The linearized equations of motion are 

3 

given by 

• • A 

^ - -20)^^ + + N^/J^ 

•• 

d) » 01 t[i + N /J (22) 

^ c y y 

0 - -30)^0 + N /J 

c z z 

where iji’, <(), and 9 refer to the yaw, roll, and pitch modes, respectively, 

OJ^ is the orbital angular rate, £2^ is the n^^ modal frequency, is 

the non-dimensional modal amplitude, and J are the principal plate 

X, y, z 

moments of inertia. 

The roll and yaw equations of motion are coupled to each other 
and the characteristic equation shows a double pole at the origin indi- 
cating instability in the roll-yaw motion. However, for an initial con- 
dition of iji(o) = 4»(o) “ 0, the roll and yaw motions will not build up. 

To study the effect of solar radiation disturbance, a square plate whose 
fundamental frequency is ten times the orbital frequency is considered. 

j ' • 

Only the first three flexible modes are included in the dynamic model 

with initial conditions of 0.01 in each mode. The transient response of 

the plate under the infliience of solar radiation pressure is shown in Fig. 4.6, 


The torque about the normal to the plate due to the first modal amplitude 
acts In one direction only (Fig. 4.3b); as the solar Incidence angle 
changes in the orbit, it is seen that the cyclic contribution due to 
averages to zero. This torque Induces a steady drift In the roll angle 
(- 1.5° In 6 orbits). The yaw motion is seen to be oscillating with 
a very small amplitude (0.3°). The solar radiation pressure disturbance 
also induces a small amplitude (0.03°) pitch oscillation. The modal oscilla- 
tions are unaffected In the presence of the solar radiation disturbance. 

The magnitude of the pitch, roll and yaw angular motions due to the 
solar radiation pressure are small becuase of the stabilizing gravity - 
gradient forces acting on the plate. 
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Fig. 4.1. A Thin Plate Exposed to Solar Radiation Pressure. 
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Fig. 4.2. The First Five Mode Shapes of a Square Plate. 
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Fig. 4.4 Force Distribution due to the Interaction of Solar Pressure 
on a Completely Absorbing Plate Vibrating in the First Mode. 



Fig. 4.5 A Square Plate in Orbit Nominally Along the Local Vertical. 











V. 


SOLAR RADIATION HEATING EFFECTS ON THIN RECTANGULAR PLATES AND BEAMS 


V.l Introduction 

The thermal gradients Induced in large flexible space structures due 
to solar radiation heating can cause large deflections depending on the 
thermal properties of the material and the geometric shape of the struc- 
ture. The thermal deflections are functions of both time and elastic 
1 2 

displacements.* The heat being radiated from a thin solid is also 
proportional to the fourth power of the surface temperatures. Therefore, 
in order to study the effects of solar heating, a number of realistic 
assumptions, depending on the problem, will have to be made to obtain 
expressions for the important thermal effects. 

The objective of the study here is to find expressions for thermal 
deflections of large, flexible, thin beams and plates exposed to solar 
heating. The major assumptions made here are: (a) the reflected solar 

radiation by the Earth (albedo) can be neglected; (b) the inherent time 
lags in the heat transfer process are very small compared with the orbi- 
tal period; (c) the radiation from the edge surfaces can be neglected; 
and, (d) the beams and plates have uniform thickness and thermal pro- 
perties resulting in a uniform temperature distribution of the surface 
facing the sun.- The effects of the Earth's shadow and shadowing due to 
a part of the structure ci*e not Included in the study. 
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V.2 Equilibrium Temperatures of Thin Plates and Beams 


The one dimensional beam bending analysis also applies to the two 

3 

dimensional plate bending. Pig. 5.1 shows the beam e.., ^ d to solar 
radiation. The solar Incidence angle, 0^, is taken to be a constant 
during a small interval of time. During this Interval the surface facing 
the sun, s^, attains temperature T^^ and the surface away from the sun, s^, 
attains a temperature, T 2 » The equilibrium temperatures, T^^ and T 2 » can 
be determined by writing the thermal balance equations. The total heat 
leaving the beam from the two surfaces, s and s., should be equal to 
the heat received by the beam. Therefore, 

(5.1) 


4 4 

E.CTT: + E„OT, «= a G cos 9, 
1 I Z Z s 1 


where , 

Ej^ and E 2 sre the emissivitles of the surfaces, and Sjj^, respectively 
a «* Stefan-Boltzman constant 

= 56.7x10^^ KW/m^ 

a *■ absorptivity of the surface, s 

d u 

2 

G » intensity of solar radiation ~ 0,8 KW/m 
The heat flowing through the plate, at equilibrium, is also equal to the 


heat radiated from the surfa ■, 

A 


E2OT2 = K(Tj^'T2 )/t^ 


(5.2) 


where , 


K •» thermal conductivity (KW/m k) of the plate material 


t " thickness of the plate 
c 


Equations (5.1) and (5.2) can be reairan^i' as 


EpOt 

Xl , ^2 -t T 2 


a G cos 0 4 . 
8 ^ 

E2CT 


E, 


4 


(5.3) 

(5.4) 
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Eqs. (5.3) and (5.4) can now be solved to obtain and T 2 by 
assuming an approximate value of or T 2 and then through niamerlcal 
Iteration. Assuming ■■ E 2 “ 0.05 and ® 0.2, the temperature 
difference, AT ■ Tj^-T 2 » Is obtained as a function of the solar Inci- 
dence angle, 6^, and various parameter ratios of shown In 

Fig. 5.2. A higher value of k^. Indicates a larger value of thermal 
conductivity and, hence, the temperature difference between the two 
surfaces becomes smaller. A plate of thickness 1cm and made of Polyamide 
(K ■ 0.25x10 ^KW/m°K) will have a maximum temperature difference of 2.3°K. 
The temperature gradient Is found to vary approximately linearly with 0^ 
from Fig. 5.2. Expressions for deflections of the plate as a function of 
the temperature gradient are developed, next. 


3 

V.3 Pure Bending of a Beam (or a Plate) : 


Fig. 5.3 shows a beam of length, %, and width, b. The temperature 

of the mid-plane of the beam is denoted by T^. The temperature of the 

surface facing the sun, s^, is then T^+(AT/2) , and the temperature of 

the surface, s», is given by T -(AT/2). According to the theory of beam 

n 

bending analyzed in Ref. 3, we have 


/ T y d A 

y 

where 

z is the transverse deflection of the beam, 

» coefficient of linear expansion 
J “ moment of inertia of the beam about the y axis 

y 

Eq. (5.5) is rewritten by evaluating the Integral 


(5.5) 


d z At 

■r—r =■ -a ~ » a constant 
dx* e t 

c 

The expression for the thermal deflection is then given by 

« at x^ 

Z n — ct — — 

e t 2 
c 


(5.6) 


(5.7) 


The thermal deflection can be mimlmized by selecting a material of 
low coefficient of expansion or by using a material of high thermal con- 
ductivity. An increase in the thickness of the plate will also Increase 
the temperature difference (Fig. 5.2) and also Increase the weight of 
the plate. Hence, the parameter, t^, should be as small as possible. 

The other in^jortant properties of materials not reflected in Eq. (5.7) 

4 

are the density and the cost of the material as shown in Table 5.1. 

For a beam of length 100m and thickness 0.01m, and made of polyamide 
(a low density and low cost material), the maximum thermal deflection is 
found to be approximately 7m. 
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If the beam is made of aluminum, the maximum deflection would be about 
2mm. Once a tolerable thermal deflection Is specified the material can 
be selected to meet the conflicting requirements of low density, high 
thermal conductivity, and low cost. In the next section the solar radia- 
tion pressure moment resulting from a thermally deflected beam (also 
applicable to a plate) Is discussed. 
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\ 


V«4 Effect of Solar Radiation Pressure on a Thermally Deflected 
Beam (or Plate) 


The moment expressions derived by Karymov (Eqns. 4.3-A.7) are used 
to obtain the moments acting on a thermally deflected beam. The equation 
for thermal deflection Is given by (5.7). The moment on a completely 
absorbing surface Is obtained (after Integration) as, 


where 


N = a c 6 £ b j 

a o o 0 


a , b , c are the direction cosines of the solar incident 
o’ o o 


(5.8) 


radiation vector 


6 " maximum deflection (from Eq. (5.7)) => z 


The maximum deflection, 5^, can be obtained as a function of 0^^^ by 
selecting a function to represent AT In Fig. 5.2, and then by using the 
function for AT In Eq. (5.7). The moment acting on a completely reflecting 
surface Is obtained through numerical Integration, as. 


N = -0.05 a c 6 £ j 

r 0 o o 

The corresponding moment expressions for a plate are obtained as 

N *»c 6 £b(bl + aj) 

a 0 o o o-’ 

" -0.05 c 5 £b (b 1 + a £j) 

N o o o o 

r 


(5.9) 


(5.10) 


The moment on a beam (or a plate) with a surface whose general coefficient 
of reflectivity is e^, can then be obtained by using Eq . (4.8). 
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Table 5.1. Properties of Representative Materials 
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Material 

Density 

(Kg/m3) 

Expansion 
Coefficient, ttg 
(m/m °C) 

Thermal 
Cond . K 
KH/m- 

= Cost 
($/Kg) 

Graphite 

1.5x10^ 

8.3x10"^ 

8.65x10"^ 

500 

Beryllium 

1.8x10^ 

3.5x10"^ 

12.25x10”^ 

10,000 

Aluminum 

2.7x10^ 

2.1x10“^ 

28.8x10"^ 

1.1 

Polyamide 

1.13x10^ 

25x10“^ 

2.45x10"^ 

15 


AT (°K) 




Fig. 5.1. Thermal Gradient in a Beam Due Lo Solar Radiation 
Heating. 



0 30 0^ 60 90*^ 


Fig. 5.2. Thermal Gradient in a Beam as a Function 
of Solar Incidence Angle and Thermal 
Conductivity, 
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Fig. 5.3. Beam Bending Due to Solar Radiation Heating 



VI. EVALUATION OF HOOP/COLUMN COUPLING COEFFICIENTS 


The generic mode equations and the equations of rotational motion of 
a flexible orbiting body contain coupling terms between the rigid and flexi- 
ble modes and terms due to the coupling within Che flexible modes Chat are 
assumed to be small and, thus, are usually neglected when a finite element 
analysis of the dynamics of the system is undertaken. In this Chapter 
a computational algorithm that permits the evaluation of the coefficients 
in these coupling terms in the equations of motion as applied to a finite 
element model of the Koop/Column system is developed. 

Using a Newton-Euler approach, one-can express the equations of motion 
of an eieoaental mass of the system, in the frame moving with the body, as^ 

{eT + r+2{u+T+wxr+wx(0Jxr) }pdv ■ {f+e+L(q)/p}pdv (6.1) 

cm 


where p 


mass per unit volume, 


e • external forces per unit mass, 

q a elastic transverse displacements of the 
element of volume. 


£ ■ force due to the gravity on the unit mass, and 

L * the linear operator which when applied to q yields 
the elastic forces acting on the element of volume 
considered. 


r • position vector of element dv 

u) “ inertial angular velocity of the body frame 


i 

f 
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VI. 1 EQUATIONS OF ROTATIONAL MOTION 

The equations of rotational motion of the body are obtained by taking 
the moments of all the external, internal and inertial forces acting on 
the body, i.e., from Eq. (6.1) 


/ r X fa + r + 2wxr + U)xr + tjjx(ti)xr) 1 pdv 
cm 

“ / rx[L(q)/p+ f + e] pdv 


( 6 . 2 ) 


one can obtain the following form for the equations of rotational 
motion. 

-L V j. r = f 

"R 


R + r + Z = G„ + Z + C 


(6.3) 


n=l 


n=l 


n=>l 


where R = /[r x(u)xr ) - (r *w) (uxr )] pdv 

Y O O O O 

Q ” /[r xq+2r^x(a)xq)+r^x(ojxq)+qx(o)xr^) 

n*i V 

-(r -w) (wxq)-(q*w) (tuxr )] pdv 
o o 

! d''‘> - pdvxU,.-T„) + ? QlX P‘1'' 

n-1 '' 


cm o' n n y o 


G„ = / r xMr pdv 
R o o 

? G “ /[r xMq+qxMr ] pdv 

I V ® ® 


C “ /rxe pdv 

V 


r - r^+q 


M - matrix operator^ which when applied to r yields gravity-gradient forces 


a ■ acceleration of the center of mass 
cm 


= force/mass due to gravity at the undeformed center of mass 
(tt)“ modal shape vector for the n'^^ mode 




frequency of the n*-^ mode 


A^ = time dependent modal amplitude function 


6.2 




VI. 2 GENERIC MODE EQUATIONS 


P^fiEfS 

^ ouAirry 


The generic mode equation Is obtained by taking the modal components 

of all Internal* external and Inertial forces acting on the body, l.e., 

•• • • 

[a +r+2wxr4tijxr+<ijx(a)xr)] pdv 

Cl (q)/p)+f+e] pviv (6.4) 

V 

The generic mode equation Is obtained In the following form: 

A -Hd^A +^/M + ? CP /M = [g + ? g +E +D']/M (6.5) 

a n n n n n n n*-l nm n n . n 


9 

where " /C<l>^°^*wxr + <t ‘wxCoixr )] 

u Y O O 

« • 

S CP ” *wxq+ $ •Ujxq4- $ •ajx(o)xq) ] pdv 

m-1 "" V 

g ■ pdv; ? ^ 't 'Mqpdv; 

n o mn y 

B " /$ ^'^^•epdv and D' = / t ^'^^pdv*(a -f ), 
n ^ o V cm o 
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VI. 3 CARTESIAN OOMPCNENTS OF -'HiE DIFFERENT COUPLING TERMS 


The expressions for R/ Q 


s ?;(n) ;; (n) ^ 




— / ^ Cartesian components are presented in this 

xon n mn 


‘9 

section. - 

One can express the following vectors in 
their Cartesian component form as 






r(n)_. (n)" ,(n)" (n) 


n“l 


- Q^^^i + 


and 




where i, j, k are unit vectors along the body principal axes 


of inertia in the undeformed state; l, $„/ are the co- 

X ^ z 


ordinates of a point in the undeformed state. 

With the use of the component forms of the 
vectors given above, one can expand the various vector ex-’ 
pressicms given in Bqp,“( 6 ~,S) ’and ’( 6 ’.S) 


fj to + (J “J )cti ti) I i 

‘ X ;« 2 y y 2 * 


* tJyiy + 5 


+ [JjUj + (Jy-Jx)“x“yi " 
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(n) „ ” ,„(n) „(n) 
X “ “zy 
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) + 2 a [ )<d - 

' n yy zz x yx y 




zx z 


n 


- '«*x 

•^8«> + Vz">x?’«yS’> + ‘“H'layS’ 


♦H <”> ) ) 


(6,7) 




( 6 . 8 ) 


* “ 3 l‘®x?’+“yx * + 2 Mj 3 (H^>-h '^>)1 .( 6 . 9 ) 

<^ =* W + (li + w 

x' yz zy ' "y' z:< xz ' z' xy yx ' 


(n) 


"*“y‘®^’*®y* > * V*‘“y”’^®^y “z“x'«zx 


- ( 0 ^ ) 

'"z' XX yy ' 

^ ^ , IT (nm) T(»n)| 4,,. /T (“n) _r tran) , 

^mn ^\'“x^^yi "^zy ^ * “y^^zx ^xz ^ 

+ u, • + A r; 

z' xy Y X yz zy 


( 6 . 10 ) 


y ' zx 


+ i (l'™’-!.*””’) 

z xy yx 
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+ a. CO 

X y . xy yx 

+ (u,<a 

2 X zx xz 


) + 0. (L<"")+L<"”’) 

y z yz zy ' 

' X yy zz ’ 


'n 


'wn 


Y zz XX ' *' XX yy 

fe "»«’ «aa 

■ \ ll 4T' ««B 


( 6 . 11 ) 


where H 


(n) 


a0 


^ V (r\\ 

ct» $ » X, Y, 2 or 1# 2/ 3, When a is x in 


or L 


(mn) 
a 6 


the corresponding value of a in M^g is 1. In a similar 
way when a is y in or , a is 2 in M^g and when 


(n) 

‘a6 


a is 2 in or , ct is 3 in M 


(mn) 
'aS 


o8* 


The same rea- 


soning holds for 3 also. 


The expressions for and are 

y z 

obtained by the cyclic permutation of x, y, z in the 
expression for in Eq. (6.7) and the expressions for 

and are obtained by the cyclic permutation of 


X/ Y* 2 expression for G 


(n) 


in Bq . 


(6-. 9). 
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For a discretized model the expressions for the volume Integrals are 


replaced by the following summations: 
k 


H 


(n) 


a3 

T (mn), 

a 3 




(a, 3 - X,y,z) 


k 

E 

1-1 


/•A (ni) V /a (n) V 

h ^*3 ^l“l 


( 6 . 12 ) 

(6.13) 


where 

k = total number of discrete masses 
1 = Index Identifying a nodal point 

m^ “ mass concentrated at the 1*"^ node. 

^ = coordinates of m^j^ In the undeformed state 
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VI. 4 EVALUATION OF COUPLING COEFFICIENTS IN THE EQUATIONS OF MOTION AS 
APPLIED TO A FINITE ELEMENT MODEL OF THE HOOP/COLUMN SYSTEM 

VI. 4.1 Model Description 

The structural dynamic modeling of the Hoop/Column antenna has gone 
through many stages before reaching the single surface model which will 
be analyzed in this chapter. 

Initially, it had 231 nodes distributed as follows: 192 nodes on 

the 8 support circles Including the hoop (24 nodes on each circle spaced 
at 15° Intervals); 28 nodes on the mast and the feed mast; and 11 nodes 
at the points of location of the solar panels (upper and lower) , the 
S band reflector, and the feed panels (up-link and down-link)-see Figs. 

5.1 and 5.2. After reduction the number of nodes was diminished to 114 
including a total of 96 nodes on the circles: 1100, 1200, 1300, and 1400; 
7 nodes on the mast and the feed mast; and 11 nodes at the locations of 
the solar panels, the S band reflector, and the feed panels (Fig. 6.1). 


VI. A. 2 Approximate Mass Distribution 

2 

From an unpublished document prepared by the Harris Corporation, 

and submitted by NASA Langley Research Center, it has been possible to 

arrive at the mass distribution shown in Table 6.1. 9803.0 lb. out of 

the total weight of the Hoop/Column Antenna (10,070 lb.) were distributed 

between the final grid points. The distribution was done in agreement with 

the information found in the Harris Corporation document. The page numbers 

appearing in Table 6.1 refer to particular mass/moment of Inertia calcula- 

2 

tions in the Harris Corporation document. 

The small (2%) discrepancy between the calculated total mass (9803.0 lb.) 
and the stated weight of the system (10,070 lb) is thought to be attributed 
to: (1) uncertainties in the weight of specific stringers; (2) uncertainties 

inhe'ent with the finite element reduction technique where the initial mass 
mus., b? redistributed between a reduced, final number of grid (node) points; 
nci Cf) other miscellaneous uncertainties, such as the exact weight/location 
of the optical Instrument, etc. 

VI. 4.3 Cartesian Coordinates of all the Nodal Points in the Final 
NASTRAN Output 

Reference 2 contains the cylindrical coordinates of all the nodal points 
on the mast, the feed mast, and at the location of the panels and electronics. 
It also contains the Z coordinates of the planes which contain the circles 
along with their respective diameters. Thus, the Cartesian coordinates of 
all the nodal points were obtained by a simple transformation from cylindrical 
to Cartesian coordinates. 
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VI. 4. 4 Development of a Computational Algorithm for Evaluation of 
the Coupling Coefficients 

After receipt of the tape containing the modal functions, this 
information was stored in our IBM 360 in such a manner that when one 
calls subroutine^GETMP(2) , he can refer to the component of the 
mode shape vector at the grid point J by VECMP(I,J,K) . Based on this, 
an algorithm described in the flow diagram. Fig. 6.2, was designee and 
tested. As indicated in Fig. 6.2, the available data, such as: the 

Cartesian coordinates of the grid points on the mast, the feed mast and 
the ones at the locations on the appendages; and such as the mass concen- 
trations at all the nodal points are input into the software routine and 
these data will consequently have to be updated according to any develop- 
ment in the Hoop/Column modeling. The subroutine, DCS, (given the radius 
of the circles and the Z component of their centers) computes the Cartesian 
coordinates of the nodal points on the circles. 

Subroutine GETMP(2), which makes the available^ is called and 

the values of components of the desired mode shape vector at the particular 
grid point are incorporated into a loop mathematically described by Eqs. 

(6.12) and (6.13). It should be noteu that, for reasons of effectiveness, 
each coefficient is evaluated separately on the circles and on the other 
grid points and then combined to yield the corresponding coupling coefficient 
for the entire Hoop/Column system. 

The algorithm has been tested for two modes (the 7^^ and the S'^h) 
successfully, but only after the evaluation of the coefficients corresponding 
to all the 13 modes will one be able to make positive conclusions. 


VI. 4. 5 Normalization of the Equations 


In order to make a valid comparison between the terms in the rigid 
modes and the ones due to the coupling between the flexible modes, the 
equations of motion and the generic modal equations are non-dimensionallzed 
in this section. 

The Euler angles and their rates are assumed to have the following 


form; 


8 ■ TT- sin fit + 6 cos J^t and 0 = 6 cos ;it -6 fi sin fit 
U o 0 o 

When 0^-0, the amplitude of 0 is governed by that of In order 

to guarantee the small amplitude approximation for these angles, 6/(2 

must have as a maximum the order of 1/10. For a gravity stabilized 

structure, is proportional to the orbital angular velocity, 

with the proportionality constant depending on the differences in 

appropriate moments of inertia, and whether pitch, roll, or yaw motion 

is under consideration. For this reason we select 6 = w /200/3 and 

c 

• • 

similarly for <1)^ and (It is well known that the frequency of a 

rod shaped dumbbell librating near the local vertical is /3 u^). 
According to an assumed Euler sequence C(l)<t), (2)0, and (3)i|;] from 
the local vertical to the principal body axes, 
oj « (0-0) )sln ^ ^ sin 0 cos 

A C 

U) = (0-~U) )co8 (Ji + sin 0 sin (}> 

y ^ 

» ij; cog 0 + (j) 

Therefore, w = (§~w )(j> - Z 0; u Z (0-oJ ); and u) ^ ijj + and the 
X c ^ c z 

corresponding gravity gradient matrix operator, is calculated as 

.2 


M ■ to 


2 0 0 

0-10 
0 0-1 
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The terms in the rigid modal equations are normalized by dividing 
2 2 

each by S. <jiy (where i " 100m is selected as the characteristic length 

of the system and is the frequency of the 7th mode); the terras in 

2 

the generic modal equations are normalized by dividing by The 

elastic displacements are assumed in such a manner that « ***n^n’ ^n - 
2 

0 ) A . Examples of normalization follow: 
n n 

. aJ(M33-M22)CH<”>+ h'">) 

+ M33(H^> + + 2M23CH^> - 

- “zf * “7“, 


+ u CO + W W + <0 CO 

X y^ xy yx ’ y z ' yz zy z x zx 


xz 


- CO^ - (0^ ^ t- H^"^}/AC0^ 

X yy zz y zz xx z xx yy / 


In connection with the normalization of cj)^ and similar terms involving 
time derivatives of the angular velocities, it should also be noted that 
d d 

dt ” ^7 "dr" ^ dimensionless time parameter. 


VI. 4.6 Comparison of the Terms 


Table 6.2 shows the comparison between the magnitudes of the components 

of R and the components of the largest coupling coefficients in Q for 

both assumed values for A Im 1mm for all modes. 

n 

It appears that when the system is operating within the mission speci- 
fications (deflections of the order of mm), the finite element assumptions 
are valid; when the deflections are of the order of meters the coupling 
between the flexible and rigid modes should be incorporated into the equations 
of motion. 

In Table 6.3, a comparison of pertinent terms in the generic modal 

equations is given. “ A^+ is compared with ^ 

13 , 

E }/M^^ where is the modal mass associated with the n^ mode. 

m-7 

The comparison, after normalization of the different terms in the equations, 
shows that the time dependent amplitude of the modes can be approximated 
as an harmonic oscillator at least up to the point where A^ = Im for all n. 

One can therefore conclude that even though the time dependent modal 
amplitude function A^ can be modelled as an harmonic oscillator, the coupling 
between rigid and flexible modes has to be taken into consideration once 
the displacements exceed the order of mm, based on the parameters of the 
Hoop/Column system. 
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Given X,Y,Z (cartesian coordinates of the nodal 
points on the mast and feed) and the mass concen'- 
tratlon at each nodal point in the system (approxi- 
mate mass distribution) 


Call Subroutine DCS which 
computes X,Y, and Z of the 
nodal points on the df.fferent 
circles 



Subroutine GETMP (2) 
makes available 

I = mode number 
J = node number ^ 

K « x,y,z component of 

Feed + Mast Appendages 

[ Circles 




“■cB - 5 “j 


For the total system 
t»(I) rom 

•aB Sb 


Fig, 6.2. Flow Diagram Describing the Algorithm Used in 
the Evaluation of the Coupling Coefficients. 










Table 6.1 Approximate Mass Distribution at Final Grid Points (Pounds) 
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Comparison of Pertinent Terms in 
the Rigid Modal Equations 




..action 

R 

n“7^n 

X 

1.48922 E-01 

2.642388 E+00 

Y 

1.49854 E-01 

-8.816196 E+00 

Z 

9.26317 E-02 

-1.480150 E+*)0 

A "I mm 
n 

X 

1.48922 E-01 

2.642388 E-03 

Y 

1.49854 E-01 

-8.816196 E-03 

Z 

9.26317 E-02 

-1.480150 E-03 


I 


i 

t 

1 


■] 



Cl ^ 
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Table 6.3 Comparison of Pertinent Terms in the 
Generic Modal Equation 


A ■ 1mm 

2 2 

Mode Number Frequency Modal Mass R htn- P /iluj- 

rad/sec ^ ^ ^ 


- 

— 

0.7489559 

153.157 

. 629344x10’^ 

1.9543x10“® 

n 

1.3692409 

5.232954 

. 210346x10“^ 

-2.10287x10"^° 

9 

1.7471481 

3.232954 

. 34248x10"^ 

8.7299x10"^ 

10 

3.2148494 ' 

0.3046446 

.115957x10'^ 

-4.901252x10"^ 

11 

4.535031 

1.992988 

.230747x10“^ 

-3.46559x10“^ 

12 

5.5926659 

723.5216 

. 350924x10"^ 

3.226x10"® 

13 

5.7942225 

0.6561203 

.376674x10"^ 

4.678578x10'*^ 


A = Im 

2 2 

Mode Number Frequency Modal Mass R /iu. P /S.u^ 

rad/sec ^ , 


10 

3.2148494 

.3046446 

.115957 

.4901252x10"^ 

11 

4.535031 

1.992988 

.230747 

. 337854x10"'^ 

12 

5.5926659 

■^23.5216 

.350924 

3.2362x10"® 
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VII. general conclusions and reccmmsndations 


The graph theory approach has been utilized to define controlla- 
bility in. terms of term rank and input-state reachability concepts and 
applied to the placement of actuators for the proposed Hoop /Column 
system. The linear quadratic Gaussian techniques appear to offer 
greater flexibility to the controls designer in attempting to meet 
performance requirements while also maintaining propellant consumption 
at desired levels. System transient characteristics are noticeably 
degraded when an actuator^ assumed to be mounted on the hoop, is then 
removed from the system. Before definitive conclusions can be arrived 
at relative to KMS pointing and surface accuracy requirements the effects 
of both sensor and plant noise and also environmental disturbances 
should be Incorporated into the existing model and such .studies have 
been initiated. 

It Is found that control laws previously designed for the case 
where environmental effects were neglected, may be inadequate to control 
the shape and orientation of very flexible beams that are exposed to 
solar radiation. For simple systems intuitive methods (suitable gain 
adjustments) of obtaining more robust closed-loop systems in the presence 
of the disturbances are Indicated, For more complex flexible systems, 
such as platforms or antenna-type structures, numerical U)G techniques 
appear to offer the most promise and such studies are continuing, 

Preliminary results evaluating the effect of thermally Induced 
deflections due to solar heating of thin structures, Indicate that 
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such deflections may give rise to appreci^-Me disturbance torques, 
depending on the thermal properties and thickness of the materials 
used. Additional research in this important area is suggested. 

A computer algorithm has been developed and used to evaluate the 
relative magnitude of coupling terms between the rigid rotational and 
flexible modes and also intra-modal coupling terms in the general 
equations of motion using the Hoop/Column mass and material properties, 
Such coupling terms are usually not Included in finite element models 
based on the Earth-based vibrational and rigid modes only. It is seen 
that when the surface deflections are of the order of mm. (vlthln 
mission specifications), that the relative magnitude of the largest 
of these coupling coefficients is at least one order of magnitude 
smaller than the principal coefficients. For surface deflections of 
the order of centimeters or even meters (still within the small ampli- 
tude deformation assumption) the coupling coefficients can be of the 
same order as the principal ones, indicating that these terms might 
have to be added to the finite element formulation if an accurate re- 
presentation of the dynamics would be required. 


